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Abstract 

The effects of two types of randomness on the behaviour of directed polymers are 
discussed in this chapter. The first part deals with the effect of randomness in medium 
so that a directed polymer feels a random external potential. The second part deals 
with the RANI model of two directed polymers with heterogeneity along the chain such 
that the interaction is random. The random medium problem is better understood 
compared to the RANI model. 

1 Directed polymers 

A long flexible elastic string, to be called a polymer, has several features of a critical system, 
showing power law behaviours without much fine tuning 012] • An important quantity for a 
polymer is its size or the spatial extent as the length N becomes large. For a translationally 
invariant system with one end (z = 0) fixed at origin, the average position at z = N is zero 
but the size is given by the rms value 

M = 0, (rW-N", (1) 

with v = 1/2, for the free case. Similar power laws can be defined in other properties also. 
In general, such exponentslike the size exponent v define the polymer universality classand 
these depend only on a few basic elements of the polymer. In addition to the geometric 
properties, the usual thermodynamic quantities, e.g. free energy (or energy at temperature 
T = 0), entropy etc., are also important, especially if one wants to study phase transitions. 

Power lawsgenerally imply absence of any length scale in the problem. Given a mi- 
croscopic Hamiltonian with its own small length scales like the bond length, interaction 
range etc, power laws occur only in the long distance limit (large N) for thermally averaged 
quantities which require summing over all possible configurations. As a result, in the long 
distance limit, these power laws become insensitive to minute details at the microscopic 
level, a feature that may be exploited to choose appropriate simplified models to describe a 
polymer. 

In thermal equilibrium, the Boltzmann distributionultimately determines the macro- 
scopic behaviour. In most cases one may avoid the issue of probability distributionbut 
instead may concentrate only on the first few, may be the first two, moments or cumulants. 
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For example, one needs to know the average energy, entropy etc and the various response 
functions which come from the width of the distribution. Thermodynamic descriptions do 
not generally go beyond that. 

In random physical systems, one encounters an extra problem. If the randomness is not 
thermalizable ("quenched"), any quantity of interest becomes realization dependent. As a 
result, an additional disorder averaginghas to be done over and above the usual thermal 
averaging for each realization. It is therefore necessary to know if and how this extra 
averaging alters the behaviour of the system. 

Critical-like systems [3] become a natural choice for studying the effect of quenched 
randomness because it is hoped that at least some of the properties would be insensitive 
to the minute details of the randomness. Since for a critical system, the influence of the 
randomness on a long scale is what matters, it transpires that the critical behaviour will 
be affected if the disorder is a relevant variable. In the renormalization group language, a 
coupling is classified as relevant, irrelevant or marginal, if, with increasing length scale, it 
grows, decays or remains invariant, because the contribution of a relevant quantity cannot 
be ignored at long distances even if it is numerically small to start with. 

With relevant disorder, there is the obvious possibility of a change in the critical prop- 
erties (e.g. new set of critical exponents). More complex situations may include emergence 
of new or extra length scales. One may recall that approach to criticality, most often, is 
described by a diverging length scale. Developing a description of the system in terms of this 
large length scale only goes by the name of scaling theory. Emergence of any new or extra 
length scale would then alter the corresponding scaling description. In case it is possible 
to change the nature of disorder from relevant to irrelevant (say by changing temperature), 
then a phase transition would occur which will have no counterpart in the pure problem. 
For the disorder-dominated phase, on a large scale, there are possibilities of rare events (see 
Appendix [XJ which necessitates a distinction between the average value and the typical 
(e.g. most probable) value. In such situations higher moments of the quantity concerned 
become important. These are some of the aspects that make disorder problems important, 
interesting and difficult. 

i (d) , (e) 

11=0 |i=l H=2 

Figure 1: (a) A random walk in d dimensions with z as the variable along the contour of the 
polymer i.e. giving the location of the monomers, (b) Directed polymer on a square lattice. 
A polymer as of (a) can be drawn in d + 1 dimensions. This is like a path of a quantum 
particle in nonrelativistic quantum mechanics, (c) A situation where both the transverse 
space (r) and z are continuous, (d) The directed polymers on a hierarchical lattice. Three 
generations are shown for 4 bonds, (e) A general motif of 2b bonds. 

The problem of a polymer in a random medium was initiated by Chakrabarti and Kertesz 
[HE] with the application of the Harris criterion. This problem has enriched our overall un- 
derstanding of polymers and random systems in general, but still a complete understanding 
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remains elusive. Not surprisingly, the look out for simpler problems that capture the ba- 
sic essence of the original complex system gained momentum. In this context, directed 
polymer splayed a very crucial role. 

Let us define the problem here. Consider a polymer where each monomer sees a different, 
independent, identically distributed random potentials. Geometrically this can be achieved 
if the monomers live in separate spaces. One way to get that is to consider the polymer 
to be a d + 1 dimensional string with the monomers in d dimensional planes but connected 
together in the extra dimension. As shown in Fig. Q this is a polymer which is directed in 
one particular direction. Hence the name directed polymer [Bl El E] • 

For a directed polymer, the size would now refer to the size in the transverse (i-directions 
and so Eq. (0) refers to the transverse size as the length in the special ^-direction increases 
For long enough chains, it is this size that matters and enters the scaling description. 

The significance of directed polymer lies in the fact that the pure system is very well 
understood and exactly solvable in all dimensions while the random problem can be attacked 
in several different ways, a luxury not affordable in most situations. 

Two types of randomness can arise in the context of directed polymers. One type would 
involve imposition of a random external potential (random medium problem). In the second 
type, the interaction (between say two chains) is random (RANI model). In the random 
medium problem, the random potential would like to have a realization dependent ground 
state which may not match with the zero-field state. In the RANI model, the randomness 
in the interaction may lead to a change in the phase transition behaviour shown by the 
polymers. These two classes are discussed separately. 

2 Outline 

We first consider the random medium problem and then the random interaction (RANI) 
model. In both cases, disorder turns out to be marginally relevant though at two different 
dimensions. The quantities of interest in a disordered system are introduced in Sec. El For 
the random medium problem, the question of relevance of disorder, the size exponent v and 
the free energy fluctuation exponent 9 are discussed in Sees. 0] and El Sec 0] deals with the 
renormalization group (RG) for the moments of the partition function and also the special 
Bethe ansatz results for the 1 + 1 dimensional problem. A Flory approach and scaling ideas 
are also summarized there. Sec. El deals with the RG approach for the free energy via the 
Kardar-Parisi- Zhang equation. The behaviour of the overlap especially near the transition to 
strong disorder phase may be found in Sec. El . We briefly mention the unzipping behaviour 
in presence of a force and summarize some recent results for the pure case in Sec. More 
on unzipping of a two chain system may be found in Sec. 19.41 These results and scaling 
arguments are then used to visualize the nature of the ground state in Sec. El Various 
aspects of the RANI model can be found in Sees. El and The question of relevance, and 
annealed versus quenched disorder in multi chain system are analyzed in Sec. El The two 
different types of randomness or heterogeneity on hierarchical lattices are studied in Sec. 
ITU1 In the last part of this section, one may find the general validity and extension of the 
Harris criterion as applied to this polymer problem. Various technical issues are delegated 
to the Appendixes. An example of rare events is given in App. El A short review of the 
pure polymer behaviour can be found in App. |BI The issue of self-averaging and some 
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recent results about it for disordered systems are surveyed in App. O The renormalization 
group approach to polymers as used in Sec. 14.2.11 is spelt out in App. [D] in the minimal 
subtraction scheme with dimensional regularization. The momentum shell RG approach for 
the free energy is dealt with in App. [EJ All the possible flow diagrams are reviewed in App. 
iFl A short introduction to the transfer matrix approach for the directed polymer problem 
is given m App. O The exact RG for the RANI model can be found in App. [H] 
On Notation and convention: 

• To avoid proliferation of symbols, we reserve the symbol f to denote an arbitrary or 
unspecified function, not necessarily same everywhere. 

• By a sample or a realization we would mean one particular arrangement of the random 
parameters over the whole system. For a thermodynamic (infinitely large) system the 
sample space is also infinite. 

• Sample averaging is denoted by [...] av while thermal averaging is denoted by (...). 

• The Boltzmann constant is set, most often, to one, = 1. 

• "Disorder" and "randomness" will be used interchangeably. 

3 Hamiltonian and Randomness 

By definition, a directed polymer is defined in D dimensions out of which one direction is 
special that represents the sense of direction of the polymer. It is then useful to show that 
explicitly by writing D = al + 1 where d is the dimension of the transverse space. One may 
consider various possible models of the underlying space as shown in Fig. ^ 

• One may consider a lattice (square lattice in the Fig. \$}p) with the polymer as a 
random walk on the lattice with a bias in the diagonal z-direction, never taking a step 
in the —z direction. The length of the polymer is then the number of steps on the 
lattice. 

• Instead of a lattice model, one may consider continuum where both the space and the 
^-direction are continuous as shown in Fig. The polymer itself may consist of small 
rods or bonds whose length provides us with a small scale cut-off. In many situations, 
it helps to consider polymers without any intrinsic small scale cutoff. 

• Quite often it is useful to consider very special lattices, e.g. hierarchical lattices as 
shown in Figs, ^1, and because of the possibility of exact analysis. Here one 
starts with a bond and then replaces the bond iteratively by a predetermined motif 
("diamond" of 2b bonds) and the process can be iterated ad infinitum. This is a 
lattice with a well-defined dimension (see Sec. lrU|) but unfortunately without any 
proper Euclidean distance. Consequently geometric properties are not useful here. 
The effective dimension of the lattice is <i c ff = (In 2b)/ In 2, if there are 2b bonds per 
motif. A directed polymer can be taken as a random walk on this lattice starting from 
say the bottom point, going up, and ending at the top end. 
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3.1 Pure case 

Taking the polymer as an elastic string, one may define a Hamiltonian 

which gives a normalized probability distributionof the position vector r at length z from 
the end at (0, 0) 

P(r,z) = ^-L^ e~<*<*> (Kd/k B T = l). (3) 

Here k B is the Boltzmann constant and T is the temperature. One can even write down the 
distribution for any two intermediate points (r,, Zj) and (r/,£/) as 



T f -r, 
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For the lattice random walk, there is no "energy" and the elastic Hamiltonian of Eq. (J2J) 
just simulates the entropic effect at non-zero temperatures. One needs to look at the lattice 
problem in case one is interested in low or zero temperature behavior. A recapitulation of a 
few properties of polymers is done in Appendix iBl 

For a polymer of length N the probability distribution gives 

(r) = 0, (r 2 ) = N, (Kd/k B T=l) (5) 

so that the transverse size of the polymer is given by 

R „ ( r 2 ) 1/2 ~ N u with v = 1/2. (6) 

The power law growth of the size of a polymer as the length increases is a reflection of the 
absence of any "length scale" in the Hamiltonian. 

3.2 Random medium 

Let us now put this polymer in a random medium. In the lattice model of Fig. ^ each site 
has an independent random energy and the total energy of the lattice polymer is the sum 
of the energies of the sites visited. In continuum, the Hamiltonian can be written as 

H = H + dz / dr r}(r, z) S(r(z) - r) = H + / dz r}(r(z), z) (7) 
Jo J Jo 

where t] (r, z) is an identical, independent Gaussian distributedrandom variable with zero 
mean and variance A > 0, 

[v(r, z)L = °> b?( r > z ) vWi z')L = A( K r - *')$( z ~ z ')- (8) 

The averaging over r\ is to be called sample averaging, denoted by [...] av (as opposed to 
thermal averaging, denoted by (...)). With this distribution of random energies, we see 
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[H] av = Hq and so the average Hamiltonian is not of much use. Disorder averaging of sample 
dependent thermal averages is to be called quenched averaging while disorder averaging done 
at the partition function level is to be called annealed averaging. 

We shall consider the situation with one end point {z = 0) fixed. Otherwise, the polymer 
may drift in the medium to locate the best possible region that would minimize its free 
energy. Such a case, eventhough formally quenched in nature, is tantamount to an annealed 
case. 



3.2.1 Partition function 

The partition function for a polymer in a random medium or potential is given by 

VIZ e~ pH . (9) 

This is a symbolic notation ("path integral" )to denote sum over all configurations and is 
better treated as a continuum limit of a well-defined lattice partition function 

Z = e~^ z) (10) 

paths 

where the sum is over all possible paths of N steps starting from r = at z = 0. It often 
helps to define the partition function such that Z{{n = 0}) = 1 to avoid problems of going to 
the continuum limit (see Eq. ©). This is done by dividing (or normalizing) Z by Z = fi N , 
Z being the partition function of the free walker with /x as the connectivity constant (= 2 
for Fig. DJ)). 

Let us define the free energy 

F = -T\nZ, (11) 

for a polymer of length N where the end point at z = N is free. A more restricted free 
energy may be defined as 

F(T,N) = -ThxZ(r,N), (12) 

when the end at z = N is at r. 



3.3 Unzipping and response 

Quite often it is useful to force a system to change its overall configuration by applying an 
external field The response function then tells us about the rigidity of the system against 
such external perturbations. E.g., a magnetic field may be applied on a magnet and the 
magnetic susceptibility is the corresponding response function. A similar applied force for 
directed polymer is an unzipping force or a pulling force applied at one end (see Fig. [2J). 
We call such a force an unzipping forcebecause of its role in unzipping of DNA-type double 
stranded polymers [TP] . There could two different ensembles. One is a fixed force ensemble 
where one applies a force at the free end z = N and studies the change in the size and shape 
of the polymer or its response. The position of the end point is given by r = —TdF(g, N) / dg. 
The second is the fixed distance ensemble where the free end is at point r and then what is 
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the force required to maintain it at that point. Using the constrained free energy, we may 
write g = — TVF(r, N). (Here we used the same notation F to denote the free energies of 
the two ensembles. The arguments and the context would distinguish the two.) The two 
ensembles behave differently in a disordered system. 1 . 





Figure 2: A directed polymer with Figure 3: A blob picture of the polymer under a 
an unzipping force. force. Though drawn as sphere, the ^-direction is 

elongated with isotropy in the transverse direction. 

If we consider the response of a directed polymer to the unzipping force, the response 
function comes from the Hamiltonian 



H = ± K Jdz(^j + Jdzr,(v(z),z)-g-J^dz 



lg 2 N 



(13) 



The disorder is Gaussian-distributed as in Eq. 
force is 



C 



T\ 



<9 2 [lnZ] f 



[{riVj) - (r i )(r j )] E 



2 Kd 

The general response function for the 

(14) 



9 = 



with i,j representing the components. It is known in statistical mechanics that the response 
of a system in equilibrium is determined by the fluctuations. 



3.3.1 Exact result on response: pure like 

By a redefinition of the variables and using the 5-correlation of the disorder in the z direction, 
we have 



[lnZ(g)] av = [lnZ(g = 0)] av + 

from which it follows that 
TN 



g*N 
Kd ' 



Cn 



Kd 



(15) 



(16) 



1 The inequivalence of the two ensembles is known also for pure case if the force is applied at some 
intermediate point 
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as one would expect in a pure system, Eq. (J3J). And there are no higher order correlations. 

Two things played important roles in getting this surprising pure-like result: (i) The 
disorder correlation has a statistical translational invariance coming from the delta function 
in the z-coordinate, and (ii) the quadratic nature of the Hamiltonian. If disorder had any 
correlation along the length of the polymer, Eq. (JT5J1 will not be valid. 

The significance of the result is that the conventional thermal fluctuation, averaged over 
randomness, superficially does not say much about the effect of disorder. We shall see later 
that this innocuous result however contains important information. 

3.3.2 Quantities of interest 

Let us list some of the quantities which are of interest for a disordered system. 

• A random system needs to be described by the probability distribution of various 
physical quantities or by the averages and moments (over realizations). The moments 
are useful, especially in absence of full information on the probability distribution and 
also for characterization of the properties of the polymer. Since there is no unique 
partition function, one of the important probability distributions would be of the 
partition function, P{Z). Any quantity of interest needs to be averaged over such a 
distribution. Similarly the probability distribution P(F) of the free energy is also of 
interest. The thermodynamic behaviour is obtained from [F] av . In case the probability 
distribution (over the realizations) of a thermal averaged quantity X becomes very 
sharp, especially in the large size limit, one may avoid doing the disorder averaging. 
This may happen for extensive quantities because of additivity over subsamples. Such 
quantities are called self-averaging. Certain aspects of self-averaging is discussed in 
Appendix O 

• The first thing to determine is the relevance of disorder. To do so, we may write 

= MizL + (z- [z] av )}i v = in [z] av + [z2] ;]~ [Z]lv + .... (17) 

This shows the importance of the variance of the partition function. If the variance 
remains small, in the limit iV — » oo, then the polymer can be described by the average 
partition function which is more or less like a pure problem. Otherwise not. We see 
that the relevance of the disorder may be inferred from the nature of the variance of 
the partition function. 

• For the partition function we may use the simple identity 

[Z n L = [e nlnZ L = exp (j2 ^n m [Q*Z)<»]$) , (18) 

where [(lnZ) m ]lv are the cumulants. In contrast to Eq. (fTTj) . it is now the fluctuations 
of the free energythat become important. One may introduce a scaling behaviour, 
namely 

[(lnZ) 2 ]W [(lnZ) 2 ] av - [\nZ]J ~ N 2 \ (19) 
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defining a new exponent Obviously, for a pure problem 9 = 0. If higher order 
fluctuations (or cumulants) do not require any new exponent, then it is fair to expect 
[(lnZ) m ]iv ~ N rn0 . This free energy fluctuation exponent is one of the new quantities 
required to describe the random system. 

A simpler form of Eq. (|18J) is the basis of the replicaapproach for disordered systems, 
namely 



;inz] £ 



lim 
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(20) 



so that to compute the average free energy we may consider a case of n-replicas of 
the original system or after averaging, an n-polymer problem with extra interactions 
induced by the disorder though an n — > limit is to be taken at the end. A few 
possible paths to take the limit for long chains are shown in Fig. 



computer /numerical (c) 
J scaling path (b) 



1/N 




Analytical approach (a) 



Figure 4: Paths for replica approach 



Nontrivial results are expected if and only if the origin in Fig. H] is a singular point 
so that the limits n — > and N — >• oo become non-interchangeable. In other words, 
the n and N dependences should be coupled so that the appropriate path is a scaling 
path like (b) in the figure 0] 

• If we demand that In [Z™]^ is proportional to N for large N, then, apart from the 
extensive term (oc nN), there will be corrections which may be assumed to involve a 
scaling variable x = nN e . For x — > oo, F{x) ~ x 1 ^ 6 so that 

In [Z n ] av = neN + an 1/e N. (N -> oo) (21) 

This is for path (a) of Fig. 0] In contrast, if we take n — > for finite N, path (c), a 
Taylor series expansion gives 

\n[Z n ] &v = nNe + an 2 N 29 + .... (n -> 0) (22) 

Eqs. 12 11221 can be used to calculate 9, the free energy fluctuation exponent. 

• So far as the geometrical properties are concerned, we first note the lack of transla- 
tional invariance for a particular realization of disorder and therefore (r) ^ 0, but on 
averaging over randomness, translational invariance will be restored statistically and 
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so [(r)] av has to be zero. One may therefore consider the size of the polymer by the 
"disorder" correlations 

Cdis = [(r 2 )] av or C dis = [| (r) |]^ with C dis ~ N 2v (23) 

We have already seen in Eq. (|16j) that the usual correlation function is the thermal 
correlation of the pure problem and has no signature of the disorder. Of course, the 
disordered averaged probability distribution P(r, N) = [Z(t, N) / J drZ(r, A r )] av is also 
of importance. 

• For the pure case (A = 0), there is no "energy", only configurational entropy of the 
polymer. But with A ^0, there may be one or more lowest energy states. The nature 
of the ground state is an important issue for any disordered problem. For the lattice 
problem, the energy is the sum of iV(— ► oo) random energies of the sites visited. But 
it is the minimization over a connected path that makes the problem difficult. If the 
low temperature behaviour of the polymer is same as that at T = 0, the fluctuation of 
the ground-state energy will also be described by the exponent 9 of Eq. Same is 
true for the size exponent also. Such a situation requires that the disorder dominated 
phase is controlled by the "zero temperature" fixed point. 

• The problem we face in a disordered system is that there is no well defined ground state 
- the ground state is sample dependent. There is therefore no predefined external field 
that will force the system into its ground state (e.g., a magnetic field in a ferromagnetic 
problem). This is a generic problem for any random system. 

But, suppose, we put in an extra fictitious (ghost) polymer and let it choose the best 
path. Now we put in the actual polymer in the same random medium but with a 
weak attraction v with the ghost. At T = 0, this polymer will then sit on top of the 
ghost. In absence of any interaction (v — > ), the polymer would go over the ghost in 
any case if there is a unique ground state, not otherwise. At non-zero temperatures 
there will be high energy excursions and how close to the ground state we are will 
be determined by the number of common points of the two polymers. This is called 
overlap which may be quantitatively defined as 

® = ^/ *(n(*)-r a (z)), (24) 

for a given sample % and then one has to average over the disorder samples, q = [<?i] av . 

In case of a repulsive interaction, the situation will be different. If there is more 
than one ground state, the two chains will occupy two different paths and there is no 
energetic incentive to collapse on top of each other when the repulsion v — > 0+. In 
such a scenario, the overlap q(v — »■ 0+) ^ q(v — > 0—). Conversely, a situation like 
this for the overlap would indicate presence of degenerate ground states. For a unique 
ground state, the second chain would follow a nearby excited path with certain amount 
of overlaps with the ground state. A little reflection shows that overlap is associated 
with the second moment of the partition function. 

Relevance of disorder 

us first see if disorder is at all relevant. 
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4.1 Annealed average: low temperature problem 

If the disorder is irrelevant, then the average partition function [Z] av (annealed average) 
is expected to give the thermodynamic behaviour. However this extra averaging of the 
partition function may lead to a violation of the laws of thermodynamics, questioning the 
validity of the annealed averaging. This happens for directed polymers. With a Gaussian 
distribution for the random energy, from Eq. (|1U[). 

[Z} av = exp(/5 2 AiV/2) exp^m/i), (25) 

so that F/N = —T(hxfi + /3 2 A/2). The entropy obtained from this partition function 

(S = —dF/dT) by definition has to be positive which requires T >T^ = yj^jj- Annealed 

averaging will not work at very low temperatures. 

This does not necessarily mean that something like a phase transition has to happen, 
because this problem will occur for any disordered system, even finite in size. However, for 
the directed polymer problem, this does signal a phase transition, though the proof comes 
from other analysis. 

4.2 Moments of Z 

The moments can be written as the partition function of an n-polymer problem with an extra 
interaction induced by the disorder. Starting from H as given by Eq. (|7J), and averaging 
over the Gaussian distribution of Eq. (|SJ), we have 

[z\ v = J viii vn 2 ...vn n e~ m \ 

Kd r N n f civ \ 2 r N 

where H n = —J dzJ2[^) ~ J dz £J 5( Yl (z) - v 3 (z)) . (26) 

This particular form can be understood in terms of two polymers, to which one may restrict 
if the interest lies only in the second cumulant of Z. These two polymers start from the 
same point and do their random walk as they take further steps. If there is a site which 
is energetically favourable, both the polymers would like to be there. The effect is like an 
attractive interaction between the two polymers - an interaction induced by the randomness. 
For the many polymer problem (for the n-th moment), there is nothing beyond two polymer 
interaction. This has to do with the nature of correlation of disorder. 

4.2.1 Bound state: two polymer problem and RG 

For the second moment, we have a two polymer problem. The analogy with quantum 
mechanics tells us that for d < 2, any binding potential can form a bound state but a 
critical strength is required for a bound state for d > 2. In the polymer language, this 
means that any small disorder will change the behaviour of the free (pure) chain for d < 
2 (disorder is always relevant) but for d > d c — 2, if /3A < (/3A) C , the chain remains 
pure-like (disorder is irrelevant). Actually in higher dimensions (d > 2) the delta- function 
potential needs to be regularized appropriately (e.g., by a "spherical" well). Such cases are 
better treated by renormalization group (RG) which also helps in making the definitions of 
relevance/irrelevance more precise. We discuss this below [T2] . 
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In short, the second term (fluctuation in partition function) in the expansion of Eq. (jl7j) 
cannot be ignored if d < 2 or if f3A is sufficiently large for d > 2. This signals a disorder 
dominated phase for all disorders in low dimensions or at low temperatures (strong disorder) 
in higher dimensions. 

4.2.2 Expansion in potential 

We do an expansion in the interaction potential and just look at the first contributing term. 
The full series can of course be treated exactly. On averaging, the first order terms drop 
out, yielding 

[Z\ v -[Z]l = J Vn, J Vn 2 e~^e-^ J dz(3 2 A5(n(z)-r 2 (z)) + .... (27) 

This is the first order term if Eq. (|26|) is used. A diagrammatic representation is often 
helpful in book-keeping as shown in Fig. Some details of the renormalization group 
approach is given in Appendix ID1 

4.2.3 Reunion 

For this two polymer problem, the interaction contributes whenever there is a meeting or 
reunion of the two polymers at a site [T^]. At the order we are considering, there is only 
one reunion but this reunion can take place anywhere along the chain and anywhere in the 
transverse direction. The second order term as shown on the right side of Fig. El can be 
thought of as two walkers starting from origin have a reunion anywhere, thereby forming a 
loop. 




Figure 5: Renormalization of the interaction. The two polymers are represented by the two 
lines of different thickness and an intersection represents an interaction. A heavy circle on 
the left hand side represents the effective interaction that is to be used for renormalization. 

The probability that two walkers starting from origin would meet at r at z is given by 
P 2 (r, z) (Eq. (jHl)) so that a reunion anywhere is given by a space integral of this probability 
which gives 

TZ Z = J drP 2 (r,z) = (4tt^)"*, with ^ = d/2. (28) 

This exponent \I> will be called the reunion exponent. The occurrence of a power law is 
again to be noted. The eventual renormalization group approach hinges on this power law 
behaviour. 

It should be noted that the value of the reunion exponent above is that of simple Gaussian 
chains. This need not be the case with interaction. For example for repulsive interaction 
between two directed polymers, \l/ = 3/2 in d — 1 but TZn ~ A r_1 (lnA^)~ 2 in d — 2. The 
Gaussian behaviour is recovered in d > 2 [T3j . 
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4.2.4 Divergences 

The contribution in Eq. (|2"7jl at the next one loop order (0(A 2 )) as shown in Fig. El involves, 
apart from some constants, an integral over the reunion behaviour given in Eq. (|28|) . This 
integral in the limit iV — > oo is 

r N 1 

/ dz ~^j2 ~ fll ~ d/2 for d > 2 but ~ N 2 for rf < 2 - ( 29 ) 

For a finite cut-off, as is usually the case, the integral is finite for d > 2, and therefore 
[(Z — [-Z] av ) 2 ] av ~ 0(j3 2 A). This however is not the case for d < 2 with d = 2 as a borderline 
case. 

The divergence we find for d < 2 comes from the large N behaviour of the probability 
distribution and is therefore ignorant of the details at the microscopic level. In other words, 
a lattice model will also show this divergence in low dimensions. This forms the basis of a 
rigorous analysis done in Ref. [Hj, but we pursue a renormalization group approach here. 



4.2.5 RG flows 

The problem of divergence we face here is ideal for a renormalization group approach. Let us 
introduce an arbitrary length scale L in the transverse direction and define a dimensionless 
"running" coupling constant 

u(L) = ((3 3 KA)L\ e = 2 - d. (30) 

The purpose behind this length scale is to choose a tunable scale at which we may probe 
the system. One may then study the RG flow of the coupling constant as the scale L is 
changed. 

In the dimensional regularization scheme we adopt here (see Appendix \D\ for details) 
the divergence seen in Eq. |2H1 are handled by analytic continuation in d. The problem 
of convergence of the integral then appears as singularities at specific d. One then tries 
to remove these divergences in e by absorbing them in the coupling constants, thereby 
renormalizing the coupling. This in a sense takes care of reunions at small scales to define 
the effective coupling on a longer scale. One then has to rescale the system to preserve the 
original length scales. With this rescaling, one ends up with a description on a coarser scale 
with small scale fluctuation effects getting absorbed in the redefined parameters. The fact 
that the process can be implemented without any need of additional parameters is linked to 
renormalizability of the Hamiltonian. 

The effect of renormalization is best expressed by the variation of the parameters or 
coupling constants with length scale. This is called a flow equation. For the problem in 
hand, the eventual flow equation 

L ^ r = (2-d)u + u\ (31) 

where the first term follows from the definition of u while the second one is from the loop. 
The magnitude of the coefficient of the u 2 term is not very crucial because at this order, 
this coefficient can be absorbed in the definition of the u itself. What matters is the sign of 
the u 2 term. General cases are discussed in Appendix |F] 



13 



d<;2 



d>2 

u=0 u=|e| 



d = 1 




u=0 



(a) 



(b) 



Figure 6: RG Fixed points for u. (a) Based on the second moment of the partition function, 
(b) Based on the KPZ equation. Arrows show the flow of u. e = 2 — d. General flow diagrams 
are discussed in App. |Fj 

Since u emanates from the variance of the disorder distribution, it cannot be negative. 
We therefore need to concentrate only on the u > part with the initial condition of 
u(L = a) = u . What we see is that for d < 2, the flow on the positive axis goes to infinity 
indicating a strong disorder phase for any amount of disorder provided we look at long 
enough length scale. An estimate of this length scale may be obtained from the nature of 
divergence for a given Uq. An integration of the flow equation gives L ~ u\j^ 2 ^ {d < 2), a 
crossover length beyond which the effect of the disorder is appreciable. 

For d > 2, there is a fixed point at u* = |e| where e = 2 — d. For u < u*, the disorder 
strength goes to zero and one recovers a "pure" -like behaviour. This is a weak disorder limit. 
But, if u > u* the disorder is relevant. 

Based on the fixed point analysis, we conclude, as already mentioned, that disorder can 
be relevant depending on the dimensions we are in (i.e. the value of d) and temperature or 
strength of disorder. In particular, one finds 

1. A disorder-dominated or strong disorder phase for all temperatures for d < 2. 

2. A disorder dominated or strong disorder phase at low temperatures for d > 2. 

For d > 2, one sees a phase transition by changing the strength of the disorder or 
equivalently temperature for a given A. This is an example of a phase transition induced by 
disorder which cannot exist in a pure case. 

It is to be noted that the phase transition (the unstable 0(|e|) fixed point) occurs because 
of the positive u 2 term in the flow equation of Eq. (j3DJ). At d = 2, u is marginal (no L 
dependence) but renormalization effects lead to an eventual growth. Such parameters are 
called marginally relevant. Any marginally relevant variable will produce an unstable fixed 
point, and hence a phase transition, in dimensions higher than the dimension in which it is 
marginal. A general statement can then be made: Disorder is expected to produce a new 
phase transition if it is marginally relevant at its critical dimension. 

The new phase transition (a critical point) is to be characterized by its own set of 
exponents. An important quantity is the length scale behaviour. The flow equation around 
the fixed point for d > 2 shows that one may define a diverging "length-scale" associated 
with the critical point as 



In the critical dimension d = 2, there are log corrections. In the weak disorder phase where 
the disorder is irrelevant, [lnZ] av ~ In [Z] av , and therefore one may put a bound on the 
transition temperature T c for a lattice model as T c > Ta as defined below Eq. (|25|). 



£ ~ \u — u 



" c , with C = 1/|2 - d\. 



(32) 
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Attempts were made to determine ( by numerical methods and verify the RG prediction. 
However, the results from specific heat ^3] and size calculations JB] agree neither with each 
other nor with the RG result of Eq. (J32)) . This remains an open problem. 

We come back to the strong disorder phase in Sec. 15.61 



4.3 Bethe ansatz and 9 

For the directed polymer problem, a mapping to a quantum problem helps in the evaluation 
of \Z n \ w at least in d = 1. For a Gaussian distributed, delta-correlated disorder, [Z"] av 
corresponds to the partition function of an n-polymer system with the Hamiltonian given 
by Eq. (|26j). Noting the similarity with the quantum Hamiltonian with z playing the role of 
imaginary time, finding iV -1 In \Z n \ w for N —> oo is equivalent to finding the ground state 
energy E of a quantum system of n particles. This problem can be solved exactly only in 
one dimension (d = 1) using the Bethe ansatz (213 121] • This gives the ground state energy 
as 

E = -K(n - n 3 ) in d = 1, (33) 
which gives 

= i (34) 

from Eq. (}2"T]) As we shall see below, this implies v = 2/3 so that the polymer has swollen 
far beyond the random walk or Gaussian behaviour. What looks surprising in this approach 
is that there is no "variance" (2nd cumulant) contribution. It is just not possible to have 
a probability distribution whose variance vanishes identically. This is a conspiracy of the 
iV — > oo limit inherent in the quantum mapping and the value of the exponent 9 that 
suppressed the second cumulant contribution (see Eq. (J2D0- 



4.3.1 Flory approach 

Using the quantum analogy, we may try to estimate the ground state energy in a simple 
minded dimensionally correct calculation based on the assumption of only one length scale. 
Such an approach generally goes by the generic name of "Flory approach". The elastic 
energy is like the kinetic energy of quantum particles which try to delocalize the polymers 
(random walk) while the attractive potential tries to keep the polymer together. For n 
polymers there are n(n— l)/2 interactions. We take the large n limit so that if the particles 
are bound in a region of size R, the energy is (using dimensionally correct form with R as 
the only length scale) 

„ n n 2 A , , 

which on minimization gives E ~ n 3 consistent with the Bethe ansatz solution. At this 
point we see the problem of the replica approach if the limit is taken too soon. Since our 
interest is eventually in n — > 0, we could have used in this argument the linear term of the 
combinatorics. That would have made energy "extensive" with respect to the number of 
particle and replaced the disorder-induced attraction by a repulsion (note the negative sign). 



15 



The end result would however have no n 3 dependence. This is a real danger and any replica 
calculation has to watch out of these pitfalls. Quite strangely we see that the correct answer 
comes by taking n — > oo first and then n — > or, probably better to say, by staying along 
the "attractive part" of the interaction only. 

4.3.2 Confinement energy 

Suppose we confine the polymer in a tube of diameter D. This is like the localization length 
argument used to justify the energy in the quantum formulation. The polymer in the random 
medium won't feel the wall until its size is comparable to that, D ~ Nq which gives the 
length at which the polymer feels the wall. Elastic energy of a blob is D 2 /Nq. But because 
of the tube, the polymer will be stretched in the tube direction. One may then consider the 
polymer as consisting of free N/Nq blobs aligning with the force, so that the energy is 



This gives the known form 1/D 2 used in the quantum analogy, Eq. ()35|) . ( and consistent 
with dimensional analysis) but for v = 2/3, this gives 1/D. 

A cross-check of this comes from the energy of a blob. Each blob has the fluctuation 
energy Nq and so free energy per unit length F/N ~ Nq/N ~ D~ 2 ( l ~ u )/ U . 

5 Analysis of free energy: specialty of directed poly- 
mer 

Another unique feature of this directed polymer problem is that there is a way to study 
the average free energy and implement RG directly for the free energy bypassing the n — > 
problem of the replica analysis completely, giving an independent way of checking the results 
of replica approach. 

5.1 Free energy and the KPZ equation 

For a polymer, the partition function satisfies a diffusion or Schrodinger-like equation. This 
equation can be transformed to an equation for the free energy F(r, z) = —T In Z(r, z). This 
is the free energy of a polymer whose end point at z is fixed at r. To maintain the distance 
fixed at r a force is required which is given by g = —S7F. If we want to increase the length 
of a polymer by one unit, we need to release the constraint at the previous layer (think of 
a lattice). The change in free energy would then depend on the force at that point, and of 
course the random energy at the new occupied site. The change dF(r)/dz being a scalar 
can then depend only on the two scalars V ■ g and g 2 . A direct derivation of the differential 
equation for the free energy shows that these are the three terms required. The differential 
equation, now known as the Kardar-Parisi-Zhang equation [7] , is 




(36) 




(VF) 2 + ri(r,z). 



(37) 
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If we can solve this exact equation and average over the random energy 77, we get all the 
results we want. 

One may also write down the equation for the force in this "fixed distance" ensemble as 
ds. T 1 

af=2x Vg -2x g ' Vg + v " (r ' 2 >- < 38) 

This equation is known as the Burgers equation. 



5.2 Free energy of extension: pure like 

We want to know the free energy cost in pulling a polymer of length N from origin (where 
the other end is fixed) to a position r. For the pure case, the free energy follows from Eq. 
(jnj) (with K inserted) as 

7 2 

F(r, N) - F(0, N) = -K T -. (pure) (39) 

The probability distribution for the free energy can be obtained by choosing g = Kdr-^/N 
in Eq. (fT3*j) as. 

P(F(T„)) = P(F(0) + -£f). (40) 
One then gets the surprising result, 

A 2 

[F(v, N) - F(0, iV)] av = g K Jf- (border) (41) 

Therefore, on the average the stretching of a chain is pure-like (elastic) with the same elastic 
constant though the fluctuation is anomalous (9 ^ 0). This is analogous to the pure-like 
result for the correlation function, Eq. (|16j). These results have a far reaching consequence 
that in a renormalization group procedure, the elastic constant must remain an invariant 
(independent of length scale). As we shall see, this invariance condition puts a constraint 
on v and 9, making only one independent. 



5.3 RG of the KPZ equation 

To analyze the nonlinear KPZ equation, an RG procedure may be adopted. This RG is 
based on treating the nonlinear term in an iterative manner by starting from the linear 
equation. This is a bit unusual because here we are not starting with a "Gaussian" polymer 
problem, rather, a formal linear equation |2S] that does not necessarily represent a polymer. 
Leaving aside such peculiarity, one may implement the coarse-graining of RG to see how the 
couplings change with length scale. 



5.3.1 Scale transformation and an important relation 

Under a scale transformation x — > bx, and F — > b e ^ v F, the randomness transforms 

like A — > b~ d ~^ v 'A. This transformation done on Eq. (j4*Tj) shows that for K to be an 
invariant (no 6-dependence) we must have, 

9 + 1 = 2u. (42) 
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This is trivially valid for the Gaussian pure polymer problem but gives a relation between 
the free energy ffuctuationand the size of the polymer. This is borne out by the intuitive 
picture we develop below. This relation gives the size exponent i/ = 2/3ind=l. 
The equation in terms of the transformed variables is then 

®L = ^b ( - 1 - 2 ^ V 2 F - A- b ^-2u + l)/u (VF) 2 + b (l-<h,-26)/u ^ z) (43) 

The 6-dependent factors can be absorbed to define new parameters, except for K. 

The temperature however gets renormalized as T — > Tb^ l ~ 2u ^ u . Its flow is described by 
the flow equation 

QT 1 

L— — = T (to leading order) (44) 

UIj v 

For v > 1/2, T(L) — > 0. The disorder dominated phase is therefore equivalent to a zero 
temperature problem. In other words, the fluctuation in the ground state energy and ground 
state configurations dominate the behaviour at low temperatures in situations with v > 1/2. 
It is this renormalization that was missing in Sec 14.2.11 

5.4 RG flow equation 

The nonlinearity (or g 2 ) contributes further to the renormalization of the temperature 
through the appropriate dimensionless variable u = (KA/T 3 )L 2 ~ d (same as in Eq. ()30|) . As 
in the previous section, the important flow equation is for this parameter u (upto constant 
factors). The renormalization of temperature acquires subtle ^-dependence that introduces 
a new element in the flow equation. Some details may be found in Appendix |EJ We quote 
the flow equation here as 

L tL^ 2 - d ^ + ^- < 45 > 

General cases of such a flow equation are discussed in Appendix |FJ 

One immediately sees a major difference with the flow equation Eq. (J31)) for d = 1. 
Because of the change of sign of the quadratic term in Eq. (J45|) . there is now a fixed point 
for d = 1 as shown in Fig. E(b). This flow equation does not behave properly in a range 
d G [1.5,2) but that is more of a problem of implementation of RG than directed polymer 
per se, and so, may be ignored here. Note also that no extra information can be obtained 
from Eq. (J45|) for d > 2 other than what we have obtained so far in Sec. EJ namely the 
existence of a critical point. 

This approach however has the advantage of getting the renormalization of temperature 
by u. Some details of this RG is given in Appendix |Ej Eq. f)44j) then gets modified to 



T dTfl-2v2-d\ rT1 

di = \~ir~ + ~ir u ) < 46 > 

If we now demand scale invariance at a fixed point of u, we can determine the exponent v. 
At d=l, the stable fixed point u* = 2 from Eq. (|45jl then gives the exact exponents at d — 1 

»-§.'-§■ (47) 
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in agreement with the Bethe ansatz results mentioned above. To get 9 the exponent relation 
Eq. ()42|) (from invariance of K) has been used. The RG results for d = 1 are expected to 
be exact. 



5.5 Critical point for d > 2 

For d > 2, the unstable fixed point is u* = 0(|e|), (d = 2 — e). This gives at the critical 
point v = 1/2 + 0(e 2 ) indicating the possibility of 

v = 1/2, and 9 = (48) 

to be exact. 

One may argue j2H] for 9 = in the following way. At T = T c , thermal fluctuation 
enables the polymer to get out of the trap set by the random potential ("ground state"). 
Just above the critical point, on a scale determined by the correlation length of the critical 
point, the random potential scaling is set by b e with the value of 9 at T c . As T — > T c +, 
the length scale diverges and therefore the relevant energy scale would also grow with the 
same exponent. However a critical point implies the energy scale to be 0(T C ) which is finite. 
These can be reconciled if 9 = at the critical point. One then gets v = 1/2. Though this 
is the same as that of a Gaussian polymer, we shall see later that the polymer has extra 
non-Gaussian features. 



5.6 Strong disorder phase for d > 2 

Unlike the strong disorder phase at d=l, the absence of any fixed point for the strong 
disorder phase for d > 2 in this approach forbids quantitative results about the phase itself. 
In addition, the behaviour of the strong disorder phase in d — 1 can be obtained by various 
methods. Unfortunately, there are few concrete results in higher dimensions d > 2. Most 
reliable values of the exponents come from various numerical approaches based on the KPZ 
equation. A recent estimate for d = 2 is ^7j 

1 jv = 1.67 ± 0.025, and 9 = 0.229 ± 0.05. (49) 

Numerical studies indicate that v decreases as d increases. A question arises about the 
existence of an upper critical dimension d = g?ucd such that v — 1/2 for d > g?ucd- For 
example, higher loop contributions in the RG of Sec. 14.2.51 show singularities at d = 4 which 
could indicate d = 4 as another critical dimension. Various analytical approaches jHJ UHl ITU] 
suggest ducD = 4, or even nonintegral dimension [20] • But numerical simulations and other 
arguments [21112210] suggest g?ucd = oo. This issue is yet to be resolved. 

The fact that the size exponent v (often called wandering exponent) is different from 1/2 
has important implications in various applications, especially for flux lines in superconduc- 
tors. For example, confinement of a flux line in presence of many other flux lines would 
lead to a steric repulsion |27j (similar to the confinement energy in Sec. I4.3.2J) and the in- 
teraction of the vortices may lead to an attractive fluctuation induced (van der Waals type) 
interaction |28j. 
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6 Overlap 



In a replica approach, occupancy of different ground states may be achieved by "replica 
symmetry breaking" (i.e. various replicas occupying various states) but the difficulty arises 
from the n — > limit. In the case of directed polymer, we have argued that the degenerate 
states occur only rarely and therefore the effect of "replica symmetry breaking" if any has 
to be very small. This is why the Bethe ansatz gave correct results without invoking replica 
symmetry breaking. 

The method to compute the overlap was developed by Mukherji [2H|. By introducing a 
repulsive potential v J dz5(ri(z) — t 2 (z)) for the two polymers in the same random medium, 

+ / dz r)(r 2 (z), z) + / dz v 5(r 12 (z)). (50) 
Jo Jo 

The free energy F(r 1 (z),r 2 (z), z,v) satisfies a modified KPZ type equation 



dz ^ \2K 1 2K 

i 

which can be studied by RG. The mutual interaction has no effect on the single chain 
behaviour but the interaction gets renormalized. The flow equation for the dimensionless 
parameter u of Eq. remains the same. The exponent relation of Eq. ()42|) also remains 
valid. The interaction gets renormalized as 

T 9v fi-e , u\ 

L *r = — "d+o «> (52) 



dL \ v 

where v is in a dimensionless form. For the pure problem (0 = 0, v = 1/2) this reduces to 
the expected flow equation of Eq. (J3Tj) for repulsive interaction (u — > —u). For overlap one 
needs only the first order term because we need v — > 0. 

The overlap can be written in a polymer-type scaling form q = iV s ^{vN~^ u ), where 
S = 9 — (j)v — 1. The above RG equation for v shows that the exponent E = at the stable 
fixed point for u of Figl^b) at d — 1. However, E < at the transition point for d > 2. 
This means that the overlap vanishes at the transition point from the strong disorder side 
as q ~ \T-T c \Wt. 

This approach to overlap can be extended to m-chain overlaps also, which show a non- 
linear dependence on m at the transition point [SUj- This suggests that eventhough the size 
exponentis v = 1/2 Gaussian like, there is more intricate structure than the pure Gaussian 
chain. Overlaps of directed polymers on trees have been considered in Ref. |31| . A case of 
cross-correlation of randomness (each polymer seeing a different noise) has been considered 
by Basu in Ref. 



7 Unzipping: pure case 

Unzipping was first considered in the context of DNA ^U]. However the same ideas play a 
role here. Let us consider a pure case of a directed polymer with one end fixed at origin 



20 



and with an attractive interaction with a line at r(z) = (instead of being in a random 
medium). The Hamiltonian for a delta- function interaction can be written as 

j pN / o \ 2 pN pN o 

very similar in form with Eq. (fTHj) except here we have an attractive interaction instead of 
a random medium. 

For the zero force case, there is a critical unbinding transition at v = t> c . For <i < 2, 
f c = 0. The pulling force would like to align the polymer in the direction of the force while 
the interaction would like to keep the polymer attached to the rod. At zero temperature the 
unzipping transitiontakes place at a force where the binding energy is compensated by the 
force term. Upto a geometric factor a, this is given by Nv = agN. At nonzero temperature, 
the entropic effects are to be taken into account, which may be done by using the quantum 
analogy. 

The problem can be mapped on to a quantum Hamiltonian, albeit non-hermitian, for a 
particle of co-ordinate r 

H q (g) = ^(p + tg) 2 + V(v), (54) 

in units of h(= ksT) = 1 and mass = 1, with p as momentum. For long chains N oo the 
free energy is the ground-state energy of this non-hermitian Hamiltonian. A phase transition 
takes place whenever the ground state is degenerate. The analysis done in Ref. fU] shows 
that if the ground state energy ( i.e. the binding energy of the polymer per unit length) is 
Eq, then the critical force is given by 

g c = 2^E~ ^\v- Vc \^- d \. (55) 

where the w-dependences of E close to v c , for general d, is used. In fact if the bound state 
has extensive entropy, then there is a possibility of a re-entrance at low temperatures (see 
Sec. 19.41 This however is not possible in this case in hand. 

8 Nature of ground states and excitations 

Powered by the quantitative estimates of the free energy fluctuation and size exponents, we 
now try to generate a physical picture. 

8.1 Rare events 

We have seen that there is a low temperature region (in lower dimensions for all T) where 
randomness results in a new phase but the response to an unzipping forceis the same as for 
the pure system. For the pure case as N — > oo the width of P(r^,N) increases. Hence 
the increase of Ct with N. With randomness, for T — > we need to look for the minimum 
energy path. Let us suppose that there is a unique ground state, i.e. E(tn) or F(tn) is 
a minimum for a particular path. This tells us that as the temperature is changed, T still 
low, the polymer explores the nearby region so that the probability distribution gains some 
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width which is determined by the thermal length. Susceptibility would be the width of the 
distribution and this is independent of N. This cannot satisfy the relation given by Eq. 
(|42|) . If we invoke the the unzipping argument, then we need to exert a force exceeding the 
critical force to take the polymer out of the bound state and so the response to a small force 
(g — > 0) would be insignificant. The situation will not be any better even on averaging over 
the random samples if every sample has a unique ground state. 

However, it may happen that most of the samples have unique ground-states but once 
in a while (rare samples) there is more than one ground state which happens to be far away 
from each other. Suppose there are such rare samples, whose probabilities decay as N~ K , 
where the paths are separated by N", then the contribution to the fluctuation from these 
samples would be j\f 2u ~ K . In case k = 2u — 1, we get back the exact result. The relation 
of Eq. (|42|) tells us k = 9. The rare events control the free energy fluctuation. From the 
unzipping point of view, the threshold in such rare cases is at zero force because a small force 
can take a polymer from one ground state to another one, gaining energy in the process. 
Following Ref. [321, one ma y ar g ue that the gain in energy from the force should be similar 
to the energy fluctuation. Assuming a scaling of the force, g ~ gN a , then gN a N u ~ N e 
which gives a = — u — 1/3 in d — 1. This argument implies that if the average stretching 
is proportional to g and to N, i.e., [ < r > ] av ~ gN, then one should get a linear plot if 
[ < r > ] SCV /N U is plotted against gN v ~ 6 . The surprising feature is the sample dependence of 
such a plot. For a directed polymer, these quantities can be obtained by a transfer matrix 
calculation which is exact for a given sample and finite N. Some details on the transfer 
matrix approach are given in Appendix O As shown in Fig. [7J one sees steps with an 
overall linear dependence. The susceptibility within a flat step is zero as seen in the plot of 
the fluctuations. 
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Figure 7: (a) Plot of < r > /N v vs gN v e for d — 1. Two different values of N and in 
different realizations of disorder, (b) Corresponding fluctuation in position. From Ref. [Bl] . 

What we see here is that though the average behaviour is the same as that of the pure 
system, the underlying phenomenon is completely different; the average thermal response 
is determined by the rare samples that have widely separated degenerate ground state and 
the probability of such states also decays as a power law of the chain length. This picture 
also shows the ensemble dependence. What we have discussed is the fixed force ensemble. 
In the fixed distance ensemble, if we keep the end point fixed at r and try to determine 
the force to maintain it at that point, then by definition, the force comes from a small 
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displacement around r. Such small displacements will never lead to the big jumps that 
ultimately contribute to the average susceptibility. This difference in behaviour in the two 
ensembles is one of the important features of quenched randomness. 

In a given sample elastic energy ~ r 2 /N ~ N 2v ~ l . The pinning energy would also grow 
with length say as N e . We see, 9 = 2u — 1 = 9. One way to say this is that the sample to 
sample fluctuation and the energy scale for a given sample are the same. 



These results can now be combined for an image of the minimum energy paths. If the 
end points at z = N are separated by r N u , the paths remain separated (each path 
exploring an independent disordered region) until they join at Az ~ r 1 ^, after which they 
follow the same path. If the end points are separated by a distance r ^> N u , the two paths 
explore independent regions and they need not meet. This picture (FigJHJ) is often alluded 
to as the river-basin network. 

8.2 Probability distribution 

For a pure polymer, the probability distribution of the end point is Gaussian but it need not 
be so for the disordered case. One way to explore the probability distribution is to study 
the response of the polymer as we take it out of its optimal or average position, e.g. by 
applying a force. In a previous section we used the fixed distance ensemble where the end 
point was kept fixed and we looked at the force g required to maintain that distance (see 
Eq. Here we consider the conjugate fixed force ensemble. 

8.2.1 Response to a force 

Let us apply a force that tries to pull the end of the polymer beyond the equilibrium value 
r ~ R . In equilibrium, the average size R or extension by the force can be expressed in a 
scaling form 



This is because for zero force one should get back the unperturbed size while the force 
term may enter only in a dimensionless form in the above equation where the quantities 
available are the size i?o and the thermal energy. For small g, linearity in g is expected. 
This requirement gives 




Figure 8: Various paths for various locations of the end point. 




(56) 



R = R ^9 (k B =l), 



(57) 
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R is not proportional to TV if v ^ 1/2 {Rq ~ N u ). The polymer acts as a spring with T/Rq 
as the effective spring constant. 

8.2.2 Scaling approach 

Let us try to develop a physical picture and the corresponding algebraic description (called 
a scaling theory). The polymer in absence of any force has some shape of characteristic size 
Rq. The force stretches it in a way that it breaks up into blobs of size £ 5 = T/g. For size 
< t; g the polymer looks like a chain without any force and these blobs, connected linearly 
by geometry, act as a "new" polymer to respond to the force by aligning along it. We now 
have two scales Rq and £ g = T/g. A dimensionless form is then 

R = R f^yS.g (58) 

Now each blob is of length N g so that £ g = Ng and there are N/N g blobs. We therefore 
expect 



R = W i 9 = ml- llv = N (|) (59) 



This gives a susceptibility x — dR/dg ~ g( 1 ~ 2u )/ v 



8.2.3 Probability Distribution 

Let us try to get the susceptibility of Eq. (J59|) in another way. Let us assume that the 
probability distribution for large R is 

P(r) ~ exp (-(r/R ) s ) . (60) 

The entropy is given by S(r) = — In P. The free energy in presence of a force which stretches 
the polymer to the tail region is given by 

F = T(r/R ) s -gr. (61) 

This on minimization gives 

rp / \ 6-1 

1 I r 



3 Rq \Ro) ^ ^ 

By equating this form with Eq. (j59j) . we get 

s = jh < 63 > 

and 

V(i-")' 



P(R) ~ exp 



r 
Rt] 



(64) 
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For v = 1/2 we do get back the Gaussian distribution. 

The above analysis, done routinely for polymers, relies on the fact that there is only 
one length scale in the problem, namely, the size of the polymer. If we are entitled to do 
the same for the disorder problem, namely only one scale, Rq ~ N", matters, then the 
blob picture goes through in toto. The chain breaks up into "blobs" and the blobs align as 
dictated by the force. Each blob is independent and the polymer inside a blob is exploring 
its environment like a directed polymer pinned at one end. The probability distribution is 
therefore given by Eq. (fMJ) which for d — 1 is 



P(r) ~ exp(-|r| 3 /iV 2 ) (d = 1). 



(65) 



If we use the relation AF = F(x, N) — F(0, N) ~ x 2 /N, then the above probability distri- 
bution can be mapped to the distribution of the free energy as 



P(F) ~ exp 



|AF| 



1/2(1-!/)' 



exp 



|AF| 3 / 2 



(d=l). 



(66) 



This has been tested numerically [B]. See, e.g. Ref. jHSj for more recent work. 



9 Random interaction - RANI model 

So far we have been considering the problem of random medium. A different situation arises 
if there is randomness in the interaction of polymers. This is the RANI model [361 137| . 
Consider the problem of two directed polymers interacting with a short range interaction 
as in Eq. (|HT|) but take v to be random. Such problems are of interest, especially in the 
context of DNA where the base sequence provides heterogeneity along the chain. In this 
DNA context, the randomness can be taken to be dependent only on the z coordinate and 
not on others like the transverse position r. It can be written as 



N 

H mt = I dz v [1 + b(z)} £(n(*) - r 2 (*)), 
'o 



(67) 



where the randomness in introduced through b(z). We take uncorrelated disorder with a 
Gaussian distribution 



[Kz)]„ = 0, [b(z 1 )b{z 2 )] m = AS{z 1 -z 2 ). 



(68) 



This would correspond to uncorrelated base sequence of a DNA problem. The full Hamil- 
tonian can be written as 



H 



N 



dz 



dri(z) 
dz 



dz 



N 



dz v [l + b{z)\ V{r 12 {z)). 



(69) 



where rj(z) is the c?-dimensional position vector of a point of chain i at a contour length z, 
and 1*12(2) = Ti{z) — r 2 (^)- Though written for a general potential V(r), we shall consider 
only short-range interaction that, when convenient, can be replaced by a <5-function. 
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9.1 Annealed case: two chain 



As expected, the average partition function (annealed averaging) would correspond to a 
pure-type problem. This however is not the case always as we see in Sec. 19.31 for more than 
two chains. 

A straightforward averaging of Z = J Dr\ Dr 2 exp(— H) using the probability distribu- 
tion of Eq. (|55|) defines an effective Hamiltonian 7i e s such that 



Dr x Dr 2 exp(-7Y cff ), 



(70) 



and it is given by 



H eS -~ 



N * 
dz ^ 



2 .2 

drAz) 



i=l 



dz 



+ v 



r N v 2 A r N 

/ dz V(v 12 (z)) - \- / dz V 2 {r l2 {z)). (71) 
Jo * Jo 



It appears from the above expression of the effective Hamiltonian that an attraction is gen- 
erated between the two chains. Now, since any short range potential under renormalization 
maps onto a 5 function potential, we can take the "minimal" effective Hamiltonian for (Z) 

as 



N 



dz 



dvtjz) 
dz 



( dr 2 (z) 
V dz 



N 

+ v I dz 5(r 12 (z)). 
'o 



(72) 



where vq is the reduced coupling constant which takes care of the attraction described earlier. 
We believe that the large length scale properties as described by Eq. (|T2*j) is same as that of 
Eq. (|71|). If necessary, we can restrict the strength of the disorder so that v , which represents 
the effective coupling between the two chains, is positive (i.e. repulsive interaction). Now 
the problem reduces to a relatively simple situation where the two chains interact with a 
pure 5-function interaction with a reduced coupling constant vq. The solution of this pure 
problem is known [T2] as discussed in the Appendix |DJ 



(a) 



(b) 





Figure 9: (a) Inter-replica interaction in the RANI model. The dotted wavy line indicates 
the "r"-type interaction between the pairs (1,2) and (3,4). (b) indicates a loop formed by 
the disorder induced interaction. 



9.2 Marginal Relevance of disorder 

When we consider the second cumulant of the partition function, we require four chains. On 
averaging over the disorder, a new "inter-replica" interaction is generated that couples the 
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original chains with the replica, 

pN 

#rep = -r / dz 6(r 12 (z))6(r M (z)), (73) 
Jo 

with fo = fgA and 3, 4 representing the two new chains, and Ty{z) = Ti — Vj. This is a special 
four chain interaction in the sense that this interaction favours a contact for chains 3 and 4 
at z if chains 1 and 2 also enjoy a contact there, though not necessarily at the same point 
in the transverse space (see Fig. ""JJ). In addition to this four chain interaction, the disorder 
also creates an effective two body attraction that changes the starting or bare interaction. 
So far as the effective two body interaction is concerned its effect on the long length scale 
is given by the RG flows of Eq. (}3~Tj) except that u can now be negative also. Assuming 
that we are at the critical point of this two body interaction, e.g. for low dimensions, at 
u* = 0, the effect of the disorder can be obtained from the RG flow of f . Defining the 
dimensionless coupling constant r = f Q L 2e where e' = 1 — d with r(L) as the dimensionless 
running coupling constant, the RG equation is given by [31"| 13*7] 

nv 

L— = 2(eV + r 2 ). (74) 

CLlj 

There are two fixed points (i) r = 0, and (ii) r* = — e'. A nontrivial fixed point becomes 
important for e' < 0, i.e., for d > l.See Appendix "HI for details. 

We see certain similarities of disorder or randomness becoming marginally relevant at 
some dimension: d c — 2 for the random medium problem, but d c = 1 for the RANI problem. 
A new fixed point emerges above this critical dimension. For the random medium problem, 
this implies the existence of a new phase and a disorder induced phase transition, but for the 
random interaction problem it defines a new type of critical behaviour. These results based 
on the exact RG analysis [213 [37] were later on also recovered from a dynamic renormalization 
group study 



9.3 Annealed case: three and four chains 

With pairwise random interaction, the RANI problem can be formulated for more than two 
chains also. In such cases, even the annealed averaging problem is not just the equivalent 
pure type problem j37|. The Hamiltonian for a four chain system is given by 

H = 2/ d2 E lT + / ^(l + 6(z))X;^W), (75) 

i<3 



where = ri(z) — Tj(z). After averaging of the partition function, using the Gaussian 
distribution of b(z), one gets the following effective Hamiltonian 



H cS = \]"dzj^(^} +v Q j\zj25{v ij {z)) 



8=1 V 7 UU i,3=l 

i<3 



Y 5(^)5^^) -2v 2 A dz 5 ( r ij) S ( r ki)- 

i,j,k JO i,j,k,l=l 

i<j<k i<j<k<l 



(76) 
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The remarkable feature of the effective Hamiltonian is that there are now two new types of 
attractive interaction one of which involves three chains, while the other one couples four 
chains together, similar to the 2-chain quenched problem. 

If we take a three chain system, the corresponding effective Hamiltonian will involve only 
the three chain term but no four chain interaction of Eq. (|76[). There is now the possibility 
of a disorder induced multicritical behavior. The four chain attractive interaction is marginal 
at d = 1 and so is the three chain interaction. The presence of these two marginal operators 
(at d — 1) however remain decoupled mainly because for directed polymers, higher order 
interactions (order = number of chains involved) do not renormalize lower order interactions. 
This has already been seen in the overlap problem for the random medium case in Sec. 
Therefore the resulting renormalization of the two new couplings are independent of each 
other, and, in fact they are the same by virtue of the nature of the interaction. Because 
of the four body interaction, we expect a disorder induced criticality as for the two chain 
quenched case, but here this happens for a real four chain system - no replica is involved. 

9.4 Unzipping 

If we consider a pure problem, there will be an unzipping transitionas described in Sec. 
above. The transition is first order in nature. In this particular case there is no thermal 
unbinding in one dimension but if we forbid crossing of the polymers, then a thermal un- 
binding is possible. On a lattice this model can be solved exactly with a phase boundary 
given by [2H] 



as shown in Fig. El One of the noteworthy features of the phase diagram is the re- 
entrance in the low temperature region [HHl EDj • At low temperatures, the unzipped chains 
are stretched by the force so that on a lattice their entropy is exponentially suppressed. 
Though one gains energy by unzipping, one loses both the binding energy and the entropy 
of the bound chain. If ln/i# is the entropy per site of the bound state in the ground-state 
then, the transition occurs at g c = v + T In \xb with a positive slope. In the case of a single 
chain interacting with a rod, the bound state has no entropy (/i# = 1) and therefore no 
re-entrance. This argument also shows that the transition is first order. 

It has been argued that the unzipping transitionfor the quenched averaged RANI case is 
second order [3T] . However for real DNA, it is not the quenched averages that matter. There 
is strong ensemble dependence and sample to sample variation. This has been exploited to 
identify point mutants by a comparison of the unzipping force in a fixed distance ensemble 
[12] -An experimental determination of the unzipping phase boundary for a real DNA has 
been reported in Ref. [13] 

10 Hierarchical lattice 

Useful information can be obtained from models amenable to exact analysis even if they look 
artificial. Real space renormalization group approach can be handled in an exact fashion 
for a class of tailor-made lattices called hierarchical lattice. Such lattices are constructed 




(77) 
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1.4 




Figure 10: Unzipping phase diagram (exact) for two mutually avoiding directed polymers 
with an attractive interaction v in 1+1 dimensions. Note the reentrance at low temperatures. 
From Ref. [M]. 

in a recursive fashion as shown in Fig. ^ The problem of a directed polymer in a random 
medium on hierarchical lattices has been considered in Ref. [J3JI33]. Here we consider the 
RANI problem on hierarchical lattices. As already noted, the effective dimensionality is 
d e & — (In 26)/ In 2 for a motif of 2b bonds. Two different situations can arise, the polymers 
interacting on shared bonds jlH] oi on contacts at vertices j!7j. Two cases are different. 
Aperiodic variation of the interaction has also been considered |48j . 



10.1 Randomness on sites 

Let us now consider interaction on the vertices. The problem is different from the bond case 
because by construction there are sites with large degrees. In other words, all sites are not 
equivalent, unlike the bond case where all bonds are on the same footing. 

Let be the partition function of a given realization of randomness and let = b L ^~ x 
the number of single-chain configurations, at the /xth generation. Here = 2 M is the length 
of a directed polymer. We define Z^{n) = [Z^] av /S'^ n , to factor out the free chain entropy. 
We call Z^n) the moments. For a given temperature, there is a critical value n = n c (jj) 
below which all moments are in their high temperature phase, in the thermodynamic limit 
of course (overbar indicating disorder average). In this limit, Z^in) approaches a fixed point 
value for n < n c , whereas, for n > n c , the moments diverge but with a finite "free energy" 
density f^{n) = (nL^)~ l In Z^. The approach to the thermodynamic limit can be written 
generically as 

Z,{n) l ' n = Z{nf/ n + B z (n)L~/ + ... (78) 

where Z(n) is the thermodynamic limit (/i — > oo) and Bz{n) is the amplitude of the finite 
size correction. 

For a given realization of disorder, the partition function can be written as (see Fig.l) 

Z, +1 = bZ^yZ\? ] + b(b-l)Sl (79) 

The first term originates from the configurations that require the two directed polymers to 
meet at a vertex. In contrast, the second term counts the "no encounter" cases. There are 
no energy costs at the two end points. The Boltzmann weight is random and for a Gaussian 
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distribution of energy, y m = y 171 . The moments of the partition function can be written as 
Z^ +1 (n) = b~ n J2 PnmZl(m) with P nm = (^) {b - l) n - m y^, (80) 

m=0 ^ ' 

and the initial condition Z (n) = 1 for all the moments because there is no interaction in 
the zeroth generation (one single bond). 

By iterating the recursion relations, the moments can be calculated exactly to arbitrary 
precision from which Bz{n) can be estimated. One finds that Bz(n) blows up as 

B z {n) ~ (n c — n)~ T as n — > n c — , (81) 

with r = 0.71 ± 0.02. This exponent is independent of temperature but depends on b. 

As mentioned, in the site version there are special sites with large degree. Whether such 
vertices lead to a Griffiths' type singularity needs to explored. Griffiths singularity in the 
context of random interactions has recently been discussed in Ref. |49j . 



10.2 Bond randomness: Problem with Harris criterion 

According to the Harris criterion, disorder is irrelevant at a critical point if the pure specific 
heat exponents < 0. We show that a more general criterion is required for the directed 
polymer problem |4*5] . 

We place two interacting polymers on a hierarchical lattice. They start from one end 
(bottom) and meet at the other end (top). See Fig. ^ There is an attractive interaction 
—v(v > 0) if a bond of the lattice is shared by two polymers. As in real space, here also the 
polymers undergo a binding unbinding transition for b > 2. 

Randomness is introduced by allowing the interaction energy to be random on each and 
every bond. The first model, model A, has independent random energy on all the 2b bonds. 
The randomness in the second model, model B, is taken only along the longitudinal direction 
so that equivalent bonds on all directed paths have identical random energy. Model B is a 
hierarchical lattice version of the continuum RANI model. 

The pure problem can be solved easily by a Real space renormalization approach where 
one needs only the renormalization of the Boltzmann factor y = exp(v/T). Let y n , Z n , and 
E n be the renormalized weight, partition function and energy at the nth generation. By 
decimating the diamonds the recursion relations are given by 

y n ^ = (y 2 n + b-l)/b, (82) 
Z n+1 = (Z 2 n + b-l)/b, (83) 

E n+1 = (84) 

The two fixed points of the quadratic recursion relation, Eq. (I82j) . are 1 and b — 1 of which 
the larger one is the unstable fixed point representing the transition point. Since y > 1 there 
is a transition at y c = b — 1 only for b > 2. The other fixed point y* = 1 corresponds to the 
high temperature limit. 

The length scale exponent^, and the specific heat exponent, a, can be obtained from 
Eqs. (JHl-dHl as 

( = ln 2 and a = 2 - C (85) 
s ln[2(6- !)/&]' V ; 
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Note that a < for b < 2 + y/2. It is clear that hyperscaling holds good with d — 1 and 
not the effective dimension d e s of the lattice. It is also gratifying to observe that the same 
v and a describe the finite size scaling form of E n . 

For the disordered case, the recursion relation for the Boltzmann weight can be written 

as 

= b-'iyMyM + V^V^ + ...V™vg*>) + Q> ~ l)*" 1 , (86) 

where yn is the Boltzmann weight in the nth generation for the disorder on the upper 
(j = 1) or lower (j = 2) part of the ith branch as in Fig. ^ To understand the effect of 
the disorder at the pure critical point we introduce a small disorder yn = y c + e« . The 
average of the disorder, [e] av , acts like the temperature as it measures the deviation from 
the pure critical point. The second moment is the measure of disorder. In principle, one 
should look at the variance of e , but at the pure critical point the variance would be the 
same as the second moment. Since our motivation is to find the flow of the disorder at the 
pure fixed point, we need only study the first two moments, starting from a finite and small 
variance. 

The crossover exponentfor the disorder is defined through the homogeneity of the singular 
part of the free energy in terms of the scaling fields \x\ (temperature) and \ii (disorder). 
Under decimation, the free energy behaves 

/sing(/il,^2) = III'* (87) 

defining <p which can be obtained from the RG recursion relations for the first two moments. 
This crossover exponent determines the relevance of disorder at the critical point and can 
be computed exactly for both models A and B. One finds a striking difference between the 
two models as 

- ln(2 ^ 3) (model A), and ^ (model B). (88) 



r ln[2(6 - l)/6] v " r ln[2(6-l)/6] 

For model A, <p is negative for all b > 2, implying irrelevance of disorder and <p 7^ a but it 
is equal to 2 — d^u, while for model B, <fi — a and not 2 — d&v. Since the randomness in 
model B is highly correlated, the Harris criterion is less expected to be valid here as opposed 
to model A. but it turns out to be so. 

In order to construct a general framework for checking the validity of the Harris crite- 
rion,we start with the Taylor expansion of the recursion relation of Eq. (jSfij) . 

e = g s (b)(e 1 +e 2 + ...e M ) + 0(e 2 l ) + ..., (M = 2b) (89) 

which defines g s (b). Simple arguments show that g s (b) determines a whose positivity is 
guaranteed if g s (b) > v^A/" -1 Now, suppose that the bonds are grouped in n sets with JVj 
bonds in the ith group such that the members of a set have the same randomness. Obviously 
^2 Ni = A/". Starting with a narrow distribution, the relevance of the disorder at the pure 
transition then requires 

g s (b) > (iVj 2 + N 2 2 + .... + NlY 1 ' 2 . (90) 
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Hence, the Harris criterion holds good if either 



g s (b) > max 



y/2 1 



or g s (b) < min 



y/2 1 

1 Vz, while for model B, n = 2 with N; 



(91) 



For model A, M = 2b and Ni — 1 Vz, while for model B, n = 2 with TVj = b. 

If the disordered models are classified by ± according to the sign of a, and I or R according 
to irrelevance or relevance of disorder, then the Harris criterionsuggests the existence of only 
two classes (+R) and (-1). On the other hand, the above inequalities allow special classes like 
(+1) and (-R) where the Harris criterion fails. Model A is in the (+1) class for b > 2 + \/2. 
Model B is in either the (+R) or (-1) class depending on b. It is possible to construct models 
that would belong to any of the four possible classes, especially (-R) 



11 Summary 

The behaviour of a directed polymer in a random medium in 1+1 dimensions seems to be 
well understood. There is a strong disordered phase at all temperatures for d < 2. For 
d > 2 renormalization group analysis shows a phase transition from a low temperature 
strong disordered phase to a weak disorder, pure-like phase. There are rare configurations 
with degenerate widely separated ground states, giving a contribution to "overlap", and 
strong sample dependent response to an unzipping force. 

The RANI model remains less understood compared to the random medium problem. 
Exact renormalization analysis establish the marginal relevance of the disorder at d = 1, 
indicating a disorder dominated unbinding transition in d > 1. Several features including a 
generalization of the Harris criterion for this criticality via relevant disorder and aspects of 
unzipping have been discussed. 

In both cases of random medium and random interaction, many issues still remain open. 
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APPENDIXES 



A Typical vs. average 

Consider a random variable x that takes two values 

X 1 = e a ^ and X 2 = e? N ,/3>l, (92a) 
with probabilities 

Pi = 1 - e~ N , and p 2 = e~ N . (92b) 

In the limit N — > oo, the average value [x] av — > e^~ x ^ N while the typical or most probable 
value is x = X\ with probability 1. On the other hand [lnx] av — > a.\f~N in the same limit, 
showing that [lnx] av is determined by the typical value of the variable while the moments 
are controlled by the rare events. Note that this peculiarity disappears if x has a smooth 
probability distribution in the sense of no special or rare events. 
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B Pure Polymers 



The universality of the "Gaussian" polymer is actually a consequence of the Central limit 
theorem. Suppose we construct a flexible polymer from bonds with independent probability 
distribution ^>(r) for a bond vector r. The end-to-end distance is given by R = so 
that the probability density of R is 

p(r) = / n dri ^ w<*e r * - r ) 

[^(k)] 7V exp(-2k-R)cik, (94) 



where ip(k) is the Fourier transform of ip( r ). For a symmetric distribution with finite vari- 
ance, ln^D(k) ~ 1 — Ak 2 /2..., which on integration over k leads to a Gaussian distribution. 
This is valid for a lattice model also. 

With the probability distribution of Eq. ©, the entropy in a fixed distance ensemble 
can be written as 



S{v) = S(0) --jp, which gives F(r) = F(0) + 3 #T" (95) 



1 r 2 „, , „,_, 1 TV 

l 

This identifies an effective spring constant for the polymer, namely 3T/2Rq. This spring 
like behaviour is purely an entropic effect. 

An important and general point is to be noted. The macroscopic quantity involves an 
"external" parameter like r which is scaled by Rq the characteristic long-length scale size of 
the polymer. That the long distance behaviour can be described by a single length scale is 
the basis of "scaling" approach to polymers. 

Another approach to scaling is to study the changes in the properties of a polymer as 
the microscopic variables are scaled. E.g., if we make a scale change, r — > br and 
the Hamiltonian of Eq. (J2J) remains invariant if v = 1/2. Under such a transformation, the 
size behaves as 

R (N)=b- 1 R (b 1 /»N), (96) 

so that by choosing b = N~ u , Rq ~ N", i.e., the size exponent becomes v — 1/2. In presence 
of interaction or disorder, v may not be obtained so simply and then renormalization group 
methods become useful. 



B.l Scaling approach in presence of a force 

If the polymer is now pulled with a force g, keeping the end at z = fixed, the polymer 
would align on the average with the force. The polymer can be thought of as consisting 
of blobs within each of which a polymer can be considered as free from the force whereas 
the blobs as unit form a still coarse grained model that shows stretching. This is shown 
schematically in Fig. Eland is used in Sec. 18.2.11 This description is called a "blob picture". 
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This picture essentially depends on the scaling idea that Rq is the only relevant scale for the 
macroscopic description. This can definitely be justified at the Gaussian level. The partition 
function with the unzipping force can be written as 



Z = / drP(r,iV)exp(/3g-r), 



(97) 



where P(r,N) is Eq. Q. The Gaussian integral can be done (keeping T explicitly in Eq. 
(J5J)) to obtain 



which is consistent with the idea of an effective spring constant of the polymer derived after 
Eq. (EH). 



A scaling approach to derive Eq. (jHEJ) would as follows. Let R g be the characteristic size 
of the polymer in presence of the force. Then, from a dimensional analysis point of view, this 
is similar to the zero force size Rq. From the nature of the force term, g is dimensionally like 
T divided by a length scale. Only lengthscale in the problem is Rq. Hence a dimensionally 
correct form is 



Note the absence of any microscopic scale in the above form. All microscopies go in R . For 
a linear law at small force we require f(x) ~ x (x — » 0) giving back Eq. (|98|) . One may 
rephrase this by saying that the force has a characteristic size £ g ~ T '/ g. 

If the polymer is confined in a tube of diameter D then the dimensionless variable is 
Ro/D. This will appear in the form for change in entropy or in "confinement energy". This 
is used in Sec. 14.3.21 

C Self-averaging 

Let us build a large system by adding blocks A, B, C, D, ... systematically so that an 
extensive quantity is a sum over its values on individual blocks. In case this averaging over 
blocks leads to a very sharp probability distribution, then no further disorder averaging is 
warranted, i.e., any large sample would show the average behaviour. A quantity with this 
property is often called self-averaging.This may not be the case if the distribution is broad 
especially in the sense discussed in App. A. A self-averaging quantity has the advantage 
that one may study one realization of a large enough system without any need of further 
disorder or sample averaging. For numerical simulations, the statistics of a non-self-averaging 
quantity cannot be improved by increasing the number of realizations. 

To be quantitative, let us choose a quantity M which is extensive meaning M = Nm 
where m is the "density" or per particle value. This is based on the additive property over 
subsystems M = J^Mj. For a random system we better write M = M(N, {Q}), with {Q} 
representing all the random variables. To recover thermodynamics, we want [M] av to be 
proportional to N for iV — > oo. Now, if it so happens that for large N 




(98) 



R g ~ Ro f(gRo/T). 



(99) 



M(N, {Q}) - Nm d , 



(100) 
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with md independent of the explicit random variables, then M is said to have the self 
averaging property. Note that no averaging has been done in Eq. ([100)1 . One way to 
guarantee this self- averaging is to have a probability distribution 

P(M/N) N -^?8(m d ). (101) 



This is equivalent to the statement that the sum over a large number of subsystems gives the 
average value, something that would be expected in case the central limit theorem (CLT) is 
applicable. This generally is the case if quantities like M for the sub-blocks are independent 
and uncorrelated random variables. 

For many critical systems CLT may not be applicable and self- averaging is not self- 
evident. A practical procedure for testing self- averaging behaviour of a quantity X is to 
study the fluctuations a 2 N = [A 2 ] av — [X]^ and then check if 

2 

Rx,n = t^W- 0, as iV > oo. (102) 

L Jav 

A quantity is not self-averaging if the corresponding R does not decay to zero. The central 
limit theorem would suggest Rx,n ~ N -1 , while a decay of Rx slower than this would be 
termed as '"weakly" self-averaging. 

We may then classify a quantity X, based on the large N behaviour, as follows: 



Rx,n ~ constant =>- non self- averaging 

~ iV _1 =>- strongly self- averaging 

~ N~ p with p < 1 =^> weakly self- averaging. 



Recent renormalization group arguments seem to suggest that if disorder is relevant 
then at the new (disorder-dominated) critical point thermodynamic quantities are not self- 
averaging j20]. The arguments leading to this extremely significant prediction of non self- 
averaging nature of critical quantities can be summarized as follows. 

Let ti = \T — T c (i, N)\/T c be a sample dependent reduced temperature where T c (i,N) 
is a pseudo-critical temperature of sample i of N sites with T c as the ensemble averaged 
critical temperature in the N — »• oo limit. In terms of this temperature, a critical quantity 
X is expected to show a sample dependent finite size scaling form 

Xt(T, N) = N p Q(UN 1/d( ) (103) 



where p characterizes the behaviour of [X] av at T c , ( being the length scale exponent. (E.g. 
p = , -f/d( when X is the magnetic susceptibility \. This is plausible because the critical 
region sets in when the size of the system is comparable to the correlation volume £ d which 
diverges as \T — T c \~^. The RG approach validates this hypothesis of Eq. ()103J) . especially 
the absence of any extra anomalous dimension in powers of N for Rx- Incidentally, this 
hypothesis, Eq. (|103|) . excludes rare events that may lead to Griffiths' singularity. Using 
this scaling form, the relative variance Rx at the critical point or in the critical region can 
be determined as 



R X ~[(5T C ) 2 ] N 2 ^, 



(104) 
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where [(£T C ) 2 ] is the sample average variance of the pseudo-critical temperature. A finite 
size scaling form for Rx can also be written down, but it is not required here. 

A random system can have several temperature scales, namely (T C (N)—T C ) and (T— T c ) ~ 
in addition to the shift in the transition temperature itself. For a system with relevant 
disorder, all these scales should behave in the same way so that typical fluctuations in the 
pseudo-critical temperature is set by the correlation volume In the finite size scaling 

limit £ d ~ N, and then 

{(5T C ) 2 } ~ N~ 2/dt: (relevant disorder) (105) 

An immediate consequence of this is that Rx approaches a constant as N — > oo indicating 
complete absence of self- averaging at the critical point in a random system. 

For a pure type critical point (irrelevant disorder) a < where a is the specific heat 
exponent (i.e. c ~ \T — T c \~ a ) of the pure system. In this case the fluctuation in T c is set 
by the size, i.e. 

l(5T c ) 2 } ~ N- 1 (irrelevant disorder) (106) 

so that, by using the hyperscaling relation 2 — a = d(, one gets < p = \a/Q < 1 where 
R x ~ N~ p . Hence all critical quantities in this case are weakly self-averaging. Moreover, it 
is the same power law involving a and (, for every critical quantity X no matter what its 
critical exponent is. 

These predictions have been verified for various types relevant and irrelevant disorders 
Exception to such non self-averaging behaviour with relevant disorder occurs 
if the T c distribution approaches a delta function for large lattices. In such a situation, one 
gets back strong self- averaging behaviour |54j . 

D Details of RG for polymers: dimensional regular- 
ization 

Some details of the renormalization group approach for polymers as done in Sec. 14.2.11 are 
given here [T2j. We consider the problem of two interacting directed polymers and study the 
second virial coefficient. The second virial coefficient is related to the two-chain partition 
function with all the ends free. Dimensional regularization is to be used here. 

For this appendix we take a simpler form of the Hamiltonian given be Eq. f!26|) as 

d r N 2 (Fir \ 2 r N 

H2 = lj Q dZ ^\d) ~ V °L dz8 ^)-*2{z)), (107) 

where vq is the bare or starting interaction strength. By introducing an arbitrary length scale 
L (may be the scale chosen to study the system or in momentum shell approach, this is the 
cutoff length), we may define the dimensionless variables uq = voL e / (2n) d ^ 2 and N = N/L 2 , 
with e — 2 — d. For long distance properties, we want L to be large. 

By definition, the second virial coefficient comes from the connected two chain partition 
function with all the ends free. An expansion in terms of the coupling constant would involve 
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polymer configurations as shown in Fig. This is incidentally identical to Eq. (|27|). Each 
line represents the probability of free polymer going from r, z to r', z, 

G(r', z'\r, z) = r -. 1 exp ( (108) 

as in Eq. (J3J). A crossing of the lines in a diagram represents an interaction at (r, z) which 
can take place anywhere requiring an integration over r and z 6 (0, N). 




Figure 11: The one loop diagram for two polymers. The labels denote the position and the 
z of the points. There are integrations over all these free coordinates. 

Since G is normalized, the spatial integration over the free end points lead to unity and 
so the dangling legs of the diagrams do not require any evaluation. One needs to do only 
the loop integrals. For example, the one loop diagram shown in Fig. ^2 corresponds to 

N 

dz' I dz I dr I dr' G 2 (r'z'\r< 



o Jo 



x J dn G(rz|ri0) J dr 2 G{rz\r 2 Q) J df x G(f x N\r'z') J df 2 G(r 2 iV|rYft09) 

which after integrations over the end coordinates r 1 ,r 2 ,f 1 ,f 2 reduces to an integral of the 
type 

n r z ' r r 

dz' / dz dr dr'G 2 (r'z'\rz) (110) 
o Jo J J 

so that integrations over the space coordinates lead to integrals of the type given by Eqs l28l 
123 Using d continuous variable, one can write 

dzz~- = — (111) 

with \I/ = d/2. This form is now valid for all d so that the loop integral may be written as, 
besides other constant factors, 

jyl-d/2 

2iVVA 2 — — — (112) 

with a simple pole at d — 2. V is the total volume. One then identifies d = 2 as the special 
dimension. 
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Using the above form with the singularity at d = 2, we may write the second virial 
coefficient as 



A 2 = 2nMVu L 



2-e 



u (l + (e/2)lnN + ...)^ 



(113) 



where a is a constant, and we used N €//2 = 1 + (e/2) In N + .... i The series has problem at 
d = 2. We try to absorb the divergences by redefining all the parameters in hand, in this 
particular case, only uq and N. It is clear the the source of divergence is the region when 
z, z' are close by. On a bigger scale it is these close by reunions which would contribute 
to the effective interaction seen. It is therefore natural that these divergences ultimately 
determine the RG flow of the coupling constant. 
Define a renormalized coupling constant as 

u = u + Diu 2 + D 2 u 2 + .... = uZ u . (114) 

Substitute in the expression for A 2 , and choose D J s to cancel the poles. Z u is called a 
multiplicative renormalization factor. 

We adopt the minimal subtraction scheme where the D's are chosen to subtract the poles 
and only the poles. The calculation is order by order and so, to one loop order, one cannot 
determine D 2 , which involves u 3 (two loop term ). 

Choosing D\ = — 1/e, we see that the divergence is removed to 0(l/e). Upto this order, 
it follows that 

u 2 

uq = u h ... (115) 

Since this absorbs all the divergences, we do not have to renormalize N. The divergence- free 
quantity is A 2 r(L,u, N) (R for renormalized; expressed in terms of u). Now, A 2 should not 
depend on L because L is put in by hand. This is ensured by demanding that 

r dA 2 

L— = 0, (116) 

where the factor of L in front has been put in for a dimensionless derivative operator. Written 
in a long form 



On 

A 2R (L,u,N) = 0, with f3(u) — L—. (117) 

oL 



Note that A 2 ^(u) = A 2 (uq), though, in general, additional renormalization factors (multi- 
plicative and/or additive) may be needed. 

By definition, (3{u) tells us how the renormalized u changes with the length scale and is 
called the beta function in RG. Some more algebra gives 

du du du du ( d\nu a \ 



dL du dL du \ du 

so that the variation of the coupling constant with scale L is given by 

fin i 

L— = p(u) = eu + u 2 /(27i). (119) 
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The factor of 2ir can be absorbed in the definition of u, as have often done. This equation, 
Eq. 11191 is called a renormalization group flow equation with the initial condition u = Uq 
for some L = L . It is analytic in e so that various dimensions can be handled with this 
equation. Initial condition may be taken as u = Uq for some L = L . In this particular case 
the beta function is exact to all orders. 

The zeros of the beta function are called fixed pointswhich can be of two types, stable or 
unstable. For this particular case, the fixed points are u* = and u* = — e. If we start with 
a very small u, the flow equation Ldu/dL m eu shows a growth of u with L if e > 0, i.e., if 
d < 2. This means u is a relevant variable at the noninteracting point. For the disordered 
system it translates to relevance of disorder at the pure fixed point. 

One gets a nontrivial fixed point at u* = — 2tt6 where /3(u) = which is an unstable fixed 
point if e < 0. This unstable fixed point for d > 2 represents the disorder induced critical 
point. 

For the disorder problem, u < is not meaningful, but for the interaction problem as 
in the RANI model, full range of u is allowed. The nontrivial fixed point for d < 2 is a 
stable fixed point and it describes the phase of repulsive polymers (fermion like). There is 
a critical binding- unbinding transition for d > 2 for pure attractive short range interaction. 
The unbinding transition is at zero potential for d < 2 

For the transition one can define a characteristic length £ so that at scales> £, the critical 
effects can be ignored. If we start with an initial value uq = u* + Ao, then one may say 
this crossover happens at some arbitrarily chosen value of Au = 1 say where An is the 
renormalized deviation from the fixed point. One may determine this by linearization as 

dAu . . , 
L— — =| e | Au, (120) 

CLIj 

which gives Am = A (L/L ) |e| . If we take A << 1 as a measure of the deviation from the 
critical point (like T — T c ), then a small deviation grows and goes over either to the stable 
f.p. for the unbound phase or to — oo for the bound phase depending on the starting initial 
sign. 

Setting L = £, we then get 

£=|A |- C , andC = A- ( 121 ) 



Since the beta function is exact, we have obtained the exact correlation length scale exponent 
for the binding-unbinding critical transition. 

At the critical dimension, d = 2, the flow equation is 

dAu , . x9 . 
L ^L = (An) ' (122) 

which gives £ ~ exp(l/A ). This exponential dependence of £ on A accounts for the 
divergence of v at e = 0. 
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E RG of the KPZ equation: momentum shell tech- 
nique 

We discuss the momentum-shell RG approach to the KPZ equation 
r)W T 1 

a7-27f v2F -27f' VF » 2 + "<'"->• < 123 > 

which is Eq. (pTTjl . More details on this equation may be found in Ref. The idea is to 
get the behaviour of the parameters of the differential equation in the long distance limit. 
The three parameters are T, K and A, of which K remains invariant so that 2v = 9 + 1, Eq. 
(I42J1 . is satisfied. Note that if the nonlinear term (i.e. the force term in the fixed distance 
ensemble) is absent, the differential equation becomes linear which can be solved exactly |23]. 
The RG scheme uses this exact solution for an iterative approach to tackle the nonlinear 
term. In contrast to the RG in the polymer approach of App. [DJ here the starting point is 
not a Gaussian polymer but a linear equation which need not represent any polymer. 
It is convenient to work in the Fourier space 

where A is an upper cutoff in /c-space (related to a real space short distance cutoff). The 
long chain, long distance limit corresponds to u, k — > 0. The KPZ equation can be written 
in the following form 

F(k,u) = G (k,u) V (k,u) + ^-G (k,uj) [ [ q.(k-q)F(q,u;-n)F(k-q,u;), (125) 



2K 

which suggests an iterative scheme. In this equation 

G? (k, u) = = 1 . , (126) 

or/(k, uo) Dk z — zoo 

for the linear equation with T> = T /2k. The effect of the nonlinear or the second term of the 
right hand side of Eq. ()125|) is to change Go to G (such that F = Grj). From this G~ l , the 
coefficient of k 2 in the k, u — > limit would give the effective temperature of the problem. 
G may be written as 

G(k,u;) = G (k,u) + ^Gl(k,u;)J^ (q - |) • (q + |) k • (q+ ^ 

xG (q -|,|-Q)Go(q+|,| + fi) G (q + |, -| - 0). (127) 
For u = and k — > 0, the above equation simplifies to 

G- X (k, 0) = G x (k, 0) (l - J-^- J dq q d -^j , (128) 



where 



u = ^L 2 -\ (129) 
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is the dimensionless coupling constant in this problem. Note that this is the same as in 
App. [DJ The power of q in the integral follows from power counting while the prefactor 
(d — 2) comes from the angular contributions of the dot products. In writing this form, 
the (i-dependent solid angle contribution has been absorbed in the definition of u. The 
nature of ultraviolet (small distance) divergence (as A — > oo) is apparent in the integral of 
Eq. (|128[) . If in the residual integration in Eq. ([128)1 . we perform a thin-shell integration 
between A(l — 51) < q < A and set A = 1, the effective temperature for the left-over long 
wavelength part is given by 

r* = r(i-«^«). (130) 

On rescaling x — > bx, we have k — > (1 — 5l)k where b = exp(Sl) ~ 1 + 81. The renormalized 
temperature is then 

f = h {i-2uy v T < w T < ( 1+5l i^a\_ (131] 



Combining Eqs. 11311 and I13UI one gets (51 = 5L/L) 

dL \ v Ad ) y J 

as quoted in Eq. (|46|). The renormalization of the disorder strength can be obtained from 
the definition 

[F*(k, u) F(k, = 2 A G(k, u) G(-k, -u). (133) 

By using Eqs. 11251 and ll28l for F and G and following the same iterative procedure as above, 
the flow equation for A can be determined as 

I g = (l=2»_ i i\^ (134) 



dL \ v 4 
Combining these two, one gets the flow equation for u 

du , 2d — 3 9 . 

i 5I = (2 - <i )" + — (135) 

as quoted in Eq. (|45|). The extra d — 2 factor in the flow equation for T leads to the stable 
fixed point at d — 1. If the (^-dependent u term in the flow of T in Eq. ()132j) is ignored, 
the resulting flow equation for u would be similar to Eq. 1)118)1 upto a constant factor for 
u. As already pointed out the coefficient of the quadratic term in the flow equation can 
always be scaled to 1 by a redefinition of u. There is a major difference between the RG 
of App. |D] and the RG done here. In App. [06 = 2 — disa small parameter and used as 
such in various expansions, though ultimately the equations remain valid for a wider range. 
Here however there is no small parameter and so the approximation cannot be controlled 
by choosing small e. 



41 



F Various flow diagrams 



We discuss the various possibilities of the flows for an RG flow equation represented by 



du 
dL 



eu + cu , 



(136) 



where u(L) is the running coupling constant at length scale L. The equation has two fixed 
points (fp), u* = 0, — e/c. For convenience, let us call the u = case as the "free" problem 
so that the u* = fp corresponds to the free case. This fixed point is to be called the trivial 
fixed point while the nonzero one is the nontrivial fp. 

marginal relevance 
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Figure 12: The various types of flow diagrams depending on e and c, the coefficients of the 
linear and the quadratic terms of the flow equation. Solid bullet (•) represents a stable fixed 
point while the star (*) represents an unstable fixed point. 

The behaviour of the coupling constant with length scale is determined by the signs of 
the two constants e and c, and the initial value u(Lq) at L = Lq. 

If e > 0, then it is a relevant variable at the free fixed point, while u is irrelevant there 
if e < 0. Special situations correspond to e = for which u is a marginal variable. The 
possibilities we need to consider are 
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1. A: e > 0,c > 0. Here the nontrivial fixed point is negative and stable. There is no 
fixed point for u > 0. A bare u < would then be equivalent to a state described by 
the nontrivial fixed point. For the random medium problem, u > and so a relevant 
u flows to large values. The resulting state cannot be described in this approach. If 
u represents the interaction between two polymers, then u < in Sec. 14.2.51 or App. 
ID1 represents a repulsive interaction while u > is for attraction. Hence one gets a 
stable nontrivial phase with a repulsive interaction in dimensions e > in region A. 
For one dimension (d = 1), one may associate this fixed point with a fermion (or hard 
core boson)-like behaviour. 

2. B: e < 0, c > 0. Change in stability of the fixed points. Here the nontrivial fixed point 
is on the positive side and the trivial or the free fixed point is stable. A system with 
negative u would behave on a long scale like a free system and so also for small values 
of u > 0. The nontrivial fixed point now represents a critical point so that a phase 
transition can be obtained by tuning u. Large values of u correspond to a different 
phase not accessible by this RG flow equation because the flow goes to infinity. 

3. C:e<0,c<0. The situation is similar to Region B except that the phase transition 
is now at a negative value of u and large negative u phase is not accessible. All positive 
u values are equivalent to the free case (asymptotically free). 

4. D: e > 0, c < 0. Here we see the free system is unstable while a positive u case is 
described by the nontrivial fixed point. 

5. e = 0. This is the c-axis, representing the marginal case. The second order leads to 
a growth of u > if c > 0. In this case, a is a marginally relevant variable while for 
c < it is a marginally irrelevant variable. There is no fixed point to describe the 
system. However, the general trend is that if for e = 0, there is a marginally relevant 
variable, then that variable leads to a critical behaviour (phase transition) for e < 0. 

The various possibilities are summarized in Fig. 

In all these cases, if the nontrivial fixed point is stable, then it represents a "critical" 
phase with characteristic exponents, while if it is unstable it represents a critical type phase 
transition with its own characteristic exponents. The reunion behaviour in Sec. 14.2.31 and 
the reunion exponents |1 3 j are examples of nontrivial exponents at a stable fixed point. 
The unstable fixed point will be associated with a diverging length scale with an exponent 
C = l/\e\ as in Eq. 

Regions A and C are related by u — > —u but others are distinct. For the polymer problem, 
as d (dimensionality) changes, the nature of flow in the partition function approach (Sec. 
14.2.51 or App. \B§ goes from region A to B. One sees a new criticality developing for higher 
d (d — 2 — e) via a marginally relevant variable. The RG flow for the RANI model (Sec. |§] 
or App. |HI) also belongs to this type. For the KPZ approach (Sec. 15.61 or App. |E]) one goes 
from region D for d = 1 to region B. 

G On Transfer Matrix 

The directed nature of the directed polymer problem makes it amenable to a transfer matrix 
approach. This is another feather in the cap of the directed polymer problem. If we know 
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the partition function at a point (x,z), we may construct the partition function for the 
z + 1th step because it is completely determined by the information available at the z step. 
The exact form of the transfer matrix would depend on the particular geometry used. 



(a) 





Figure 13: Lines show the possible steps from step z to z+ 1 on a squarre lattice. Backward 
steps are not allowed for directed polymers. (1) Lattice oriented by 45 degrees. The x and 
z axes are along the diagonals of the square lattice, (b) Axes are along the two directions of 
the square lattice but the polymer can take steps along the diagonals of the unit cell also. 

If we consider the geometry shown in Fig. QJd, where the x and z axes are along the 
diagonals of the square lattice, then the partition function satisfies (see Fig. IT5k ) 



For the standard geometry (Fig. 113b ) with x and z axes along the principal directions of 
the lattice but if we allow diagonal steps, then the partition function is given by 



where an extra 7 factor has been introduced to provide appropriate elasticity to the polymer. 
One may set 7 = 1 for a fully flexible polymer. The initial condition is 



where 5 a ^ is the Kronecker delta. The partition function for a chain of length N starting from 
(0, 0) can be obtained by iterating this equation. For a given realization of the randomness 
77, the partition function (and therefore any physical quantity) can be calculated exactly for 
finite N. For quenched averaging one has to average over various realizations, and this is 
where the exactness of the approach gets lost. 

As the length N of the polymer increases, the span of x also increases linearly so that 
for N — ► 00 one has to study an infinitely large matrix. This allows the possibility of phase 
transitions in a seemingly one dimensional problem. For numerical analysis, special care 
needs to be taken to keep track of the rapid growth of the partition function as the length 
increases. 

For T — > 0, the problem reduces to determination of the ground state energy. One may 
take the limit E(x,N) = —]im.T->oThiZ(x,N), but a direct approach is also possible. For 
geometry of Fig. ITSk . the energy can be obtained from 



Z(x, z + 1) = e -^ z+1)/T [Z(x -l,z) + Z(x + 1, z)\. 



(137) 



Z(x, z + l) = e -^ z+1)/T [Z{x, z) + <yZ(x -l,z) + ^Z{x + 1, z)\ 



(138) 



Z(x,0) = 5 xfi , 



(139) 



E(x, z + 1) = mm(E(x — 1, z) + r](x, z + 1), E(x + 1, z) + r){x, z + 1)), 



(140) 



so that the globally minimum energy path is 



E(N) 



mmE(x, N). 



(141) 
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Though we are considering the square lattice (in 1 + 1 dimensions), generalization to other 
lattices and higher dimensions are straightforward. Similarly, one may consider cases with 
random energies on the bonds instead of vertices. 

For the overlap problem of Sec. |Hl the partition function would satisfy in 1 + 1 dimensions 

Z( Xl ,X 2 , Z+l) = e -Mxi,z+l)+v(x 2 ,z+l))/T x 

J2 J2 Z ^ + P> X2 + ?> z ~) tt 1 ~ ki.**) + e- v/T 6 XuX2 ],(U2) 

p=±l q=±l 

where the last term involving v is the Boltzmann factor for the interaction on contact 
{x\ = If v is taken as a random quantity chosen from a predetermined distribution, 
then the same transfer matrix can be used to treat the RANI model also. 



H RG for the RANI model 

We show that the flow equation Eq. (J74)) is exact in the minimal subtraction scheme using 
dimensional regularization. More details may be found in Ref . |3Hl 137] . 
The Hamiltonian needed for [Z 2 ]1 V is 

H = \£ ^E(^r) -roJ o N dz6(r 12 (z))S(r u (z)), (143) 



where we have set the effective two body interaction to zero. For the cumulant, one need 
only consider the "connected" partition function for this Hamiltonian. It is advantageous 
to consider the Laplace transform of the ^-dependent partition function as 



oo 

sN[ry2lC 



Z= \ dN e- sN [Z 2 ]t v , (144) 
Jo 

the Laplace conjugate variable being s. 

(a) (b) (c) 



Figure 14: (a) The contributing ladder diagrams for \Z 2 \ C W . A thick line corresponds to 
a pair of chains. A wiggly line stands for an f factor in the evaluation of the diagrams. 
There are one and two loops in diagrams (b) and (c) respectively. Divergences arise from 
loop integrations. 

An expansion in the coupling constant fo can be arranged like a ladder ( "time-ordered" 
diagrams) as shown in Fig. El The individual pairs of chains are represented by thick 
lines. The horizontal wiggly lines in these diagrams stand for fo. Such a representation 
is possible because the 5 function in H2, Eq. (|73[). forces the members of a pair to have 
the same r, z coordinates. Each chain is described by the free distribution ( "propagator" ) 
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G(rf — r iy Zf — Zj) = [2n(zf — Zi)Y d l 2 exp[— (rj — rj) 2 /2(zj — z,)] with end points (17, zf) and 
(Ti,Zi). In conformity with the current usage we use the word "propagator" for the lines. 
Two chains are therefore described by 

G 2 (r, z) = (4irz)- d/2 G(r, z/2). (145) 

This G 2 is the propagator for the thick lines. At each wiggly line, connecting four chains (all 
four having the same chain length z), there are two integrations over the spatial coordinates 
of the two separate pairs of chains (thick lines). The loops formed out of the wiggly lines 
are only responsible for the divergence at d = 1. 

In order to trace the algebraic origin of the singularity, note that, by very nature of the 
interaction, the spatial integrations associated with the two thick lines are independent of 
each other. Each section of the thick lines, with z 1; z 2 as the end points, in a loop formed with 
the wiggly lines, contributes [z\— z 2 )~ d / 2 from the identity in Eq. ( |145|) . Since the interaction 
demands same z for the two thick lines, the z integrals involve [z\ — z 2 )~ d type factors whose 
Laplace transform would contribute T(l—d) with pole at d = 1. The two independent spatial 
coordinates which are left out after the successive use of the normalization J dr G(r, z) = 1, 
lead to a V 2 factor (total volume) for each diagram. The convolution nature of the z integrals, 
thanks to the time ordering, leads to a simple product of the individual Laplace transforms 
of the integrands, resulting in a geometric series for Z. 

As an example we consider the two loop diagram of Fig. IT4T c). After integration over 
the free ends, we are left with the following 

pN rzi pZ2 p 

f 3 / dzj dz 2 dz 3 G 2 (r 12 , zi 2 ) G 2 (r 23 , z 23 ) G 2 (r' 12 , z 12 ) G 2 (r' 23n z 23 ). (146) 

Jo J0 Jo J{r, r'} 

Here the subscripts denote the successive points along the thick lines while the two sets of 
chains are distinguished by the prime. In the Laplace space this becomes 

f 3 V 2 r 2 (e') (4ny 2d s^ 2+2e '\ (e' = 1 - d). (147) 

This can be generalized to arbitrary orders since only ladder type diagrams are involved. 

Defining the dimensionless coupling constant tq through an arbitrary length scale L as 
ro = foL 2e (4:7r)~ d , e = 1 — d, we write the series for Z to all orders in r$ as 



Z\, 0=0 = {An) d V 2 s- 2 L- 2 ' 



00 

n=l 



(148) 



It is clear from the above expression that there is a divergence at d — 1 at each order (> 1). 

A renormalization through minimal subtraction would require absorption of the poles in 
e' through 

r = r(l + air + a 2 r 2 + ...). (149) 

with a n = Y2p=i a n,p e ~ v an d t as the renormalized coupling constant. In such a scheme, 
a n,p{p 7^ n) terms are required to take care of the sub leading divergences. 

The geometric series of Eq. ()148j) guarantees that the removal of the leading poles is 
sufficient to remove the subleading ones. All the divergences can be absorbed by the choice 
a p = (— e')~ p which can be obtained by an explicit order by order calculation. 
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The (3 function is therefore exact to all orders in perturbation series and is given by 

P(r)=L— = 2(e'r + r 2 ). (150) 

There are two fixed points: (i) r = and (ii) r* = —e'. The bare coupling constant tq which 
originates from v^A, where A, the variance of the distribution, is strictly positive, requires 
a positive r. Therefore, the nontrivial fixed point for d < 1 in negative r is unphysical. It 
however moves to the physical domain for d > 1. 
Exactly at d = 1, e = 0, r grows with length L as 



r(L) = r(0) 



1 + 2r(0) ln^ 



r(0) being the coupling at length L Q . Hence, the disorder is marginally relevant. For d > 1, 
there exists an unstable nontrivial fixed point at r =| e | which separates two distinct regimes 
of disorder. If we start with a strong enough disorder, on the right side of the fixed point, 
it increases with length scale, going beyond the perturbative regime. This is the strong 
disorder phase. On the other hand, the left side of the fixed point is the weak disorder 
regime, since r flows to zero (the stable fixed point). 

References 

[1] P. de Gennes, Scaling concepts in polymer physics (Cornell U. Press, 1979) 

[2] M. Doi and S. F. Edwards, Theory of Polymer dynamics, ( Clarendon Press, 1988). 

[3] S. M. Bhattacharjee, Critical Phenomena: An Introduction from a modern perspective 
in "Field theoretic methods in condensed matter physics" (IOP Publishing, Bristol, UK 
2003/ Hindusthan Publishing 2001 (for India); |cond-mat/00110lTl 

[4] B. K. Chakrabarti and J. Kertesz, Z. Phys. B 44 , (1981) 221. 

[5] K. Barat and B. K. Chakrabarti, Phys. Rep. 258 (1995) 377. 

[6] T. Halpin-Healy and Y. C. Zhang, Phys. Rep. 254, (1995) 215 . 

[7] M. Kardar, G. Parisi and Y. Zhang, Phys. Rev. Lett. 56 (1986) 889. 

[8] D. A. Huse, C. L. Henley and D. S. Fisher, Phys. Rev. Lett. 55 (1985) 2924. 

[9] D. S. Fisher and D. A. Huse, Phys. Rev. B43 (1991) 10 728. 

[10] S. M. Bhattacharjee, J. Phys. A 33 (2000) L423. ( |cond-mat/9912297D 

[11] R. Kapri, S. M. Bhattacharjee and F. Seno, Phys. Rev. Lett. (2004) 248102. 

[12] J. J. Rajasekaran and S. M. Bhattacharjee, J. Phys. A24 (1991) L371. 

[13] S. Mukherji and S. M. Bhattacharjee, J. Phys. A26, L1139 (1993); Phys. Rev. E48 
(1993) 3427. 



47 



14] J. Z. Imbrie and T. Spencer, J. Stat. Phys. 52 (1988) 609. 

15] B. Derrida and O. Golinelli, Phys. Rev. A 41 (1990) 4160. 

16] J. M. Kim et al., Phys. Rev. A 44 (1991) R4782. 

17] L. Giada, A. Giacometti and M. Rossi, Phys. Rev. E 65 (2002) 036134. 

18] M. Lassig and H. Kinzelbach, Phys. Rev. Lett. 78 (1997) 903. 

19] J. K. Bhattacharjee, J. Phys. A, 31 (1998) L93. 

20] T. Blum and A. J. McKane, Phys. Rev. E 52 (1995) 4741. 

21] E. Marinari, A. Pagnani, G. Parisi and Z. Racz, Phys. Rev. E 65 (2002) 026136. 

22] C. Castellano, M. Marsili, L. Pietronero,Phys. Rev. Lett. 80 (1998) 3527. 

23] For Bethe ansatz, see e.g. J. Ziman, Models of disorder, Cambridge University Press, 
Cambridge, 1979 

24] M. Kardar, Nucl. Phys. B290 (1987) 582. 

25] S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. London A 381 (1982) 17. 

26] C. A. Doty and J. M. Kosterlitz, Phys. Rev. Letts. 69 (1992) 1979. 

27] T. Nattermann and R. Lipowsky, Phys. Rev. Lett. 61 (1988) 2508. 

28] S. Mukherji and T. Nattermann, Phys. Rev. Lett. 79 (1997) 139. 

29] S. Mukherji, Phys. Rev. E50 (1994) R2407. 

30] S. Mukherji and S. M. Bhattacharjee, Phys. Rev. B53 (1996) R6002. 

31] B. Derrida and H. Spohn, J. Stat. Phys. 51 (1988) 817. 

32] A. Basu, Phys. Rev. E62 (2000) 4675. 

33] M. Mezard, J. Phys. France 51 (1990) 1831. 

34] Rajeev Kapri, unpublished report. (2003). 

35] M. Praehofer and H. Spohn, J. Stat. Phys. 115 (2004) 255. 

36] S. M. Bhattacharjee and S. Mukherji, Phys. Rev. Letts. 70 (1993) 49. 

37] S. Mukherji and S. M. Bhattacharjee, Phys. Rev. E48 (1993) 3483. 

38] H. Kallabis and M. Lassig, Phys. Rev. Lett. 75 (1995) 1578. 

39] D. Marenduzzo, S. M. Bhattacharjee, A. Maritan, E. Orlandini and F. Seno Phys. Rev. 
Lett. 88 (2002) 028102. 

[40] D. Marenduzzo, A. Trovato, and A. Maritan, Phys. Rev. E 64 (2001) 031901. 

48 



[41] D. K. Lubensky and D. R. Nelson, Phys. Rev. Lett. 85, (2000) 1572. 

[42] S. M. Bhattacharjee and D. Marenduzzo, J. Phys. A 35 (2002) L349. 

[43] C.Danilowicz et al. Phys. Rev. Letts. 93 (2004) 078101. 

[44] J. Cook and B. Derrida, J. Stat. Phys. 57 (1989) 89. 

[45] M. Kardar and E. Medina, J. Stat. Phys. 71 (1993) 967. 

[46] S. Mukherji and S. M. Bhattacharjee, Phys. Rev. E52 (1995) 1930. 

[47] S. M. Bhattacharjee and S. Mukherji, Phys. Rev. E52 (1995) 4860. 

[48] T. A. S. Haddad, R. F. S. Andrade and S. R. Salinas, J. Phys. A: 37 (2004) 1499. 

[49] Y. Kafri and D. Mukamel. Phys. Rev. Lett. 91 (2003) 055502. 

[50] A. Aharony and A. B. Harris, Phys. Rev. Lett. 77 (1996) 3700. 

[51] S. Wiseman and E. Domany, Phys. Rev. E 52 (1995) 3469; Phys. Rev. E 58 (1998) 
2938. 

[52] O. Dillmann, W. Janke and K. Binder, J. stat. phys., 92 (1998) 57. 

[53] M. I. Marques and J. A. Gonzalo, Phys. Rev. E 60 (1999) 2394; 65 (2003) 057104. 

[54] S. Roy and S. M. Bhattacharjee, |cond-mat/0409012| 

[55] A.-L. Barabasi and H. E. Stanley, Fractal concepts in surface growth, (Cambridge U 
Press, UK,1995). 



49 



